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We study a class of three dimensional exactly solvable models of topological matter first put for- 
ward by Walker and Wang [arXiv:1104.2632v2]. While these are not models of interacting fermions, 
they may well capture the topological behavior of some strongly correlated systems. In this work 
we give a full pedagogical treatment of a special simple case of these models, which we call the 3D 
semion model: We calculate its ground state degeneracies for a variety of boundary conditions, and 
classify its low-lying excitations. While point defects in the bulk are confined in pairs connected 
by energetic strings, the surface excitations are more interesting: the model has deconfined point 
defects pinned to the boundary of the lattice, and these exhibit semionic braiding statistics. The 
surface physics is reminiscent of a v = 1/2 bosonic fractional quantum Hall effect in its topological 
limit, and these considerations help motivate an effective field theoretic description for the lattice 
models as variants of bF theories. Our special example of the 3D semion model captures much of 
the behavior of more general 'confined Walker- Wang models'. We contrast the 3D semion model 
with the closely related 3D version of the toric code (a lattice gauge theory) which has deconfined 
point excitations in the bulk and we discuss how more general models may have some confined and 
some deconfined excitations. Having seen that there exist lattice models whose surfaces have the 
same topological order as a bosonic fractional quantum Hall effect on a confining bulk, we construct 
a lattice model whose surface has similar topological order to a fermionic quantum hall effect. We 
find that in these models a fermion is always deconfined in the three dimensional bulk. 



I. INTRODUCTION 

It has been known for decades that two is a privileged 
dimension for point particles. In greater than two dimen- 
sions, point particles are restricted to be either bosons or 
fermions; but for D = 2, a vast number of other any- 
onic statistics are possible - some of which are believed 
to be realized by fractional quantum Hall states^. The 
study of the 2D topological phases of matter that possess 
such anyonic excitations has grown into a substantial en- 
terprise in condensed matter physics, spurred on by the 
ultimate goal of creating topologically protected qubits 2 . 

Our substantial understanding of topological phases in 
two dimensions naturally leads to the question of whether 
any of their defining features can also be realized in three 
dimensional systems. Though statistics between point 
particles in D — 3 must be fermionic or bosonic, interest- 
ing statistical interactions can still exist between vortex 
defects, between point particles and vortex defects (as 
in discrete gauge theories), or between point-like defects 
with additional structure^'. Thus one class of 3D topo- 
logical phases contains particles that have unconventional 
statistics, as they are not strictly point-like. A second 
category of 3D phases that share some of the character- 
istics of 2D topological phases are fractional topological 
insulators^]. When time reversal symmetry is present, 
the surfaces of these systems have protected gapless exci- 
tations. If time reversal is broken, however, the surfaces 
will exhibit a fractional Hall effect. These gapped sur- 
faces might thus be expected to support chiral anyonic 
excitations. 

In the present work, we focus on a third class of 3D 
topological phases, which can be realized by a family of 
Hamiltonians introduced by Walker and WangP. Specif- 



ically, we discuss phases realized by a subset of these 
models, which we will call confined Walker- Wang mod- 
els. These are reminiscent of fractional topological insu- 
lators in that the 2D surfaces of these 3D systems dis- 
play the physics of a chiral anyon model (whose topo- 
logical properties are those of a fractional quantum Hall 
state). They differ, however, in that a fractional topo- 
logical insulator must preserve time-reversal symmetry in 
the bulk- which requires a topological ground-state de- 
generacy in periodic boundary conditions^. The confined 
Walker- Wang models have non-degenerate ground states 
on a system without boundaries, and explicitly broken 
time reversal symmetry in the bulk. This breaking of 
time-reversal ensures that unlike fractional topological 
insulators^, our models do not have protected gapless 
boundary modes. Rather, they describe a system resem- 
bling a confined phase of the fractional topological in- 
sulator, in which time-reversal invariance has been lost, 
and the surface modes have been gapped such that they 
are always found in their fractional Hall state. 

The exactly solvable Walker- Wang Hamiltonians we 
study here are 3D cousins of the 2D string-net models 
introduced by Levin and WerP, and operate on the same 
underlying principles. These families of models are not 
intended to describe any known physical system; rather, 
they are interesting because they can be solved exactly, 
and capture the long- wavelength physics of certain topo- 
logical states of matter. They thus provide a framework 
in which to study the physical properties of these sys- 
tems, including their ground state degeneracy and the 
properties of their low-lying excitations. 

Our key findings are as follows. As mentioned above, 
and unlike both discrete gauge theories and true frac- 
tional topological insulators, the confined Walker- Wang 
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models that we consider have a unique ground state on 
any closed 3D lattice (i.e. a lattice that has no bound- 
ary). Thus in the absence of boundaries, the models are 
not topologically ordered in the usual sense. The possi- 
ble excitations in the bulk are line defects with an energy 
cost per unit length: in other words, the bulk admits no 
deconfined point particles. Unlike vortex lines in 3D sys- 
tems, which must form closed loops, these line defects 
can end at points in the bulk. The name "confined" 
Walker- Wang model stems from the fact that these point 
particles are confined — that is they must sit at the end- 
points of line defects that have a fixed energy cost per 
unit length (analogous to flux-tubes in a confining gauge 
theory). If the lattice does have a boundary, however, 
the system has the more conventional hallmarks of topo- 
logical order: with appropriate boundary conditions, the 
surfaces admit deconfined anyonic excitations, and the 
ground state is degenerate. Though much of our work 
focuses on the simplest (abelian) example of the family 
of confined Walker- Wang models, we will show that these 
properties also apply to models that have deconfined non- 
abelian anyons at their surfaces. 

Readers well-versed in the quantum Hall effect should 
note that surface states of the confined Walker- Wang 
models always have the topological properties of a 
bosonic quantum Hall state. However, Walker-Wang 
models with fermionic surface states also exist; in ad- 
dition to the surface topological order, these have a de- 
confined fermion and Zi topological order in the bulk, 
which we will also discuss in some detail. 

Unlike topological insulators, whose single surface 
Dirac cone cannot be realized in any 2D model with the 
same symmetries, the surface states of confined Walker- 
Wang models are no different in their long-wavelength 
properties from phases that can be realized by purely 
2D systems. Indeed, by adding an appropriate 2D layer 
to the boundary of a confined Walker- Wang model, it is 
possible to destroy the topological order resulting from 
these surface states. By this metric, they do not repre- 
sent genuine 3D topological phases of matter. 

The appearance of a chiral anyon model at the sur- 
face of an exactly solvable 3D Hamiltonian is nonethe- 
less striking. Simulating a general chiral Chern-Simons 
theory locally on a lattice is long-outstanding challenge. 
That is, though a chiral Chern-Simons theory can emerge 
as the effective description of a purely 2D system, there 
is no "fixed-point" lattice model (in the sense of Rcf. 
[9j) describing the limit in which the correlation length 
is less than the lattice spacing. Indeed, a chiral Chern- 
Simons theory has topologically protected gapless bound- 
ary modes, meaning that when the 2D system has a 
boundary, the correlation length at the boundary must 
be infinite. Thus it is believed^ that such fixed-point 
Hamiltonians do not exist for chiral 2D Chern-Simons 
theories (or, more generally, chiral 2D anyon models). 
Confined Walker- Wang models, whose Hamiltonians are 
exactly solvable, realize chiral anyon models as surface 
states of fixed-point Hamiltonians (whose correlation 



length is 0). The possibility of constructing a theory 
in terms of the boundary of a higher dimensional model 
is analogous to the domain- wall solution of the fermion 
doubling problem^!] 

The outline of the remainder of this work is as follows. 
The bulk of the paper will focus on the simplest topolog- 
ical models - those that can be described as loop gases. 
In Sec. |ll]we review the 2D loop gases: the toric code^ 
and the 2D doubled semion (DSem) model in the form 
introduced by Levin and WerP. (More expert readers 
may wish to skip this introductory section). The gener- 
alization of these quantum loop gas models to three di- 
mensions is the content of much of this paper. In Sec. |III| 
we present the 3D generalizations of these loop gas mod- 
els: the 3D Toric code^ (essentially a 3D Z2 lattice gauge 
theory) and a 3D analogue of the semion model (3DScm), 
which is our paradigm for confined Walker- Wang models. 
After describing these Hamiltonians, we discuss their low- 
energy physics for periodic boundary conditions (i.e., if 
the system is on a three-torus, T ), showing that whereas 
the toric code has deconfined point particles and a de- 
generate ground state, the 3DSem model has a unique 
ground state and no deconfined point-like excitations, 
just like a confined fractional topological insulator. 

In Sec. [IV] we study the low-energy physics of these 
models with open boundaries. We find that for 3DSem, 
boundaries introduce new ground state sectors, and that 
deconfined anyonic excitations can exist on the bound- 
ary. Sec. [V] discusses the field theory that captures the 
properties of the ground states and (confined) excitations 
of 3DSem. 

A large part of this paper focuses on the 3D quantum 
loop gases (the 3D Toric code, and the 3DSem model). 
However, these constitute only a small subset of the pos- 
sible Walker- Wang models. In Sec. |VI| we discuss how 
our results extend to the more general case. We show 
that 3DSem is one member of a large family of confined 
Walker- Wang models for which all point defects are con- 
fined in the bulk, and the ground state degeneracy is com- 
pletely determined by the topology of the boundary of the 
system. Most generally, the surface excitations of con- 
fined Walker- Wang models may be non-abelian anyons. 
We discuss the effective field theory for these lattice mod- 



els in Sec. VI D In Sec. VII we study one family of the 
even more general case: Walker- Wang models that de- 
scribe fermionic systems, whose surface states have the 
topological order of fermionic quantum Hall systems and 
whose bulk has a deconfined fermionic point excitation. 

A number of technical aspects of the models are ad- 
dressed in the appendices. Some readers will find Ap- 
pendix [C] in which we give exact expressions for the 
Hamiltonian and quasi-particle creation operators in the 
doubled semion model in the presence of arbitrary de- 
fects, of particular interest. 

Since the Hamiltonians we discuss were introduced by 
Walker and Wang 8; , it is worth emphasizing what is new 
in the present work. Our main contribution is to treat 
in detail the excitations of these Hamiltonians, particu- 
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larly on their two dimensional surfaces. This allows us 
to describe the full spectrum of "topological" features 
of the systems in question, and make more explicit the 
connection to field theories suggested in Ref. [51 We also 
discuss the possible boundary conditions of these mod- 
els, and under what circumstances the topological order 
associated with the surface can be destroyed. 



II. 2D LATTICE MODELS 

While 3D lattice models are the main focus of this pa- 
per, it is useful to first review the topological order found 
in their 2D analogues — starting with some very simple 
examples. Perhaps the simplest exactly solvable topo- 
logical model is the 2D toric code: a lattice model where 
each bond takes one of two possible states (which we can 
call spin-down and spin- up). It turns out that a very 
slight modification of the 2D toric code HamiltoniarP2l 
gives a different 2D topological pha se of matter known 
as the doubled semion (DSem) model 9 ^. Both the toric 
code and the doubled semion model can be thought of as 
quantum loop gases in that the bonds with down spins 
must form closed loops. While in this aspect the two 
models are quite similar, the resulting low energy theories 
are quite different: The toric code is the simplest example 
of a discrete gauge theory, whereas the doubled semion 
model is the simplest example of a doubled Chern-Simons 
theory^ — that is, a theory whose low energy spectrum 
is two identical anyon models with opposite chiralities. 
(In the doubled semion model, the anyon theory contains 
one species of anyon, which has "semionic" statistics in 
which exchanging two particles induces a complex phase 
±i to th e wav e function, whence the name). 

In Sec. II A we introduce the well-known toric cod e! 9 * 12 !, 



and then calculate its ground state degeneracies and clas- 
sify its low-lying excitations; experts may wish to skip 
this section. Then in Sec. Ill Bl we do the same for the 
doubled semion modeP, paying particular attention to 
the types of defects that arise and how they differ from 
those found in the toric code. The notation and tech- 
niques introduced in this section provide a firm founda- 
tion for our study of the 3D models in Sec. |III| 




FIG. 1. (Color online) This figure shows a region of the lat- 
tice for the 2D models, where we have represented the spin 
degree of freedom on each edge with a black dot. Bold black 
edges represent <r z = — 1 spin configurations, while dashed 
edges represent <r z = +1. The figure also indicates the edges 
involved in the definition of B v ; the three green edges in s(v) 
are acted on with a z matrices as shown. The edges involved 
in B p differ between the toric code and the DSem model. In 
the toric code B p acts on the six red edges in dp with a x 
matrices as shown, and does not act on the six blue edges in 
s(p) i.e. /(cr z ) = 1. In the DSem model, in addition to acting 
on the red edges with a x , the plaquette operator also acts on 
the six blue edges in s(p) with / (a z ) = i' 1-0 " ^ 2 . 



vertex operator B v takes value ±1 depending on whether 
there are an even/odd number of down-spins on the edges 
coming into vertex v. The plaquette operator B p flips the 
spins on every edge of a plaquette p and commutes with 
B v because it flips a pair of spins at each vertex v. Since 
clearly [B V ,B V >] = [B P ,B P '] — 0, the Hamiltonian is a 
sum of commuting operators, allowing us to write down 
its ground states exactly. 



Ground States 



A. The toric code 

The Hilbert space of the toric code consists of a two- 
state system <r 2 = ±1 on each edge of a honeycomb lat- 
tice as shown in Fig. [T] The Hamiltonian is given by 

n °?-En"r. w 

v i€s(v) P i<£dp 

B v B p 

where ^ v runs over all vertices, and ^ p over all plaque- 
ttes, s(v) is the set of three edges attached to v and dp is 
the set of six edges bounding plaquette p (see Fig. [IJ . The 



Using the fact that all B v and B p have eigenvalues ±1 
we see that a lower bound on the ground state energy is 
obtained when 

Yl af = l for all v J]_ of = 1 for all p . (2) 

Any state which satisfies all of these conditions is au- 
tomatically a ground state because it saturates a lower 
bound on the Hamiltonian. Remarkably such states ex- 
ist. Henceforth we work with a basis of states for which 
the spins on the edges of the lattice are a z eigenstates, 
and we choose a 'loop gas' representation of these states 
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FIG. 2. This figure shows the graphical rules for the toric 
code and DSem models. Row (a) represents the fact that the 
ground state only involves vertices with B v = 1. The dia- 
grams in (b)-(d) serve two purposes. Firstly, they tell us the 
relative amplitudes of loop gas configurations in the ground 
state e.g. row (c) tells us that configurations related by re- 
moving a closed loop occur with the same amplitude in the 
toric code ground state, but with a relative minus sign in the 
DSem ground state. Second, these diagrams provide a neat 
graphical mnemonic for the definitions of string operators. 



as diagrams with only down-spin edges colored. The first 
condition in Eq. Q ensures that there are an even num- 
ber of down-spins on the edges attached directly to any 
vertex, so the ground state must be made of a superpo- 
sition of spin configurations which all involve only non- 
intersecting closed loops. The second condition in Eq. ^ 
tells us that any two spin configurations that are related 
to each other by flipping all the edges of a plaquette ap- 
pear with the same coefficient in the ground state. This 
has three important consequences: Given two spin config- 
urations obeying the first constraint in Eq. Q , if the loop 
gas pictures of these spin configurations can be related 
by deforming loops over plaquettes (as shown in Figure 
2i>), by removing contractible loops (as shown in Figure 
2:), or by fusing loops together (as shown in Figure [2ji) , 
then the spin configurations appear in the ground state 
with the same amplitude. This is because the actions 
of deforming loops, removing contractible loops and fu- 
sion can all be implemented by applying a sequence of 
Bp operators, and B p = +1 on the ground state. As a 
concrete example, if four configurations differ only in a 
small region in the manner shown in Fig. [3|a)-(d) then 
they all appear in the ground state with the same ampli- 
tude. Loop deformation relates Fig. 3|a) and (b), loop 
collapsing relates Fig. ^a) and (c) and loop fusion relates 
Fig.ga) and (d). 

Now that we have concrete relations (Fig.[2jb)-(d)) be- 




FIG. 3. This diagram shows a specific small region of the lat- 
tice for four distinct spin configurations. The configurations 
are taken to be the same everywhere outside this region. For 
the toric code, all of these configurations occur in the ground 
state wave function with the same amplitude. This is because 
(b), (c), and (d) can all be made to look like (a) using loop 
deformation, loop collapsing and fusion respectively. In the 
DSem model, (a) and (b) occur with the same phase because 
they are related by deforming a loop, (a) and (c) occur with 
a relative minus sign because they differ by the presence of a 
single closed loop, (a) and (d) also occur with a relative minus 
sign because they are related by a single fusion (Fig. [2^d)). 



tween the coefficients of different spin configurations in 
the ground state we can deduce the ground state degen- 
eracy for periodic boundary conditions in both directions 
(i.e. on the torus T 2 ). Starting with any configuration of 
closed loops one can show that it appears in the ground 
state with the same amplitude as one of the four 'canon- 
ical' configurations labelled (a)-(d) in Fig. |4j These are 
distinguished by the winding number parities 
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(3) 



where c x and c y are cuts along the x and y directions 
respectively. Each distinct ground state is an equal su- 
perposition of all configurations related to a canonical 
ket by the equivalences in Fig. [2](b)-(d) ; some of the kets 
in each superposition are shown in Fig. [4] 



2. Excitations 

The Hamiltonian Eq. Q has two types of excitations: 
vertex defects where B v = —1, and plaquette defects 
where B p = — 1. These defects are exact eigenstates 
because both B p and B v are conserved (that is, they 
commute with the Hamiltonian). Two further features 
of these excitations will be important in what follows. 
Firstly, an excitation necessarily involves a pair of vertex 
or plaquette defects. Secondly, both types of excitations 
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FIG. 4. This figure shows the 4 independent ground state 
sectors in the toric code with periodic boundary conditions in 
both directions (i.e., on a torus). The 4 "canonical" kets can 
be taken to be the four spin configurations on the far left side 
labeled (a)-(d). Each for the 4 ground states is a superposition 
of all spin configurations which can be made to look like the 
canonical ket by using the equivalence rules shown in Fig. [2]d- 
d. For the toric code the relative coefficient is always +. The 
ground state structure for the DSem model is similar except 
the kets marked with a ± appear with negative coefficients, a 
fact that can easily be checked using the rules in Fig. [2] 

are 'deconfined', meaning the energy of a pair of vertex 
(or plaquette) defects is independent of their relative sep- 
aration. 

We can understand both of these features by consid- 
ering how defects are created. To create a vertex defect 
at v, one must flip an edge spin in s(v). This, however, 
creates a violation on the adjacent vertex V2 sharing this 
edge. We can remove this violation by flipping a second 
edge spin in sfa)- However, this just moves the violation 
to a third vertex V3, and so on. Hence we can never cre- 
ate a lone vertex defect but we can create a pair of vertex 
defects by flipping the spins along a continuous path C of 
edges by acting with the operator Wy{C) = Yl il£C erf, as 
shown in Fig. [5j We call this a vertex-type string opera- 
tor. Crucially the string operator commutes with all pla- 
quette operators along its length (because it only involves 
a x, s), so only has an energy cost from the vertex defects 
at its endpoints. This tells us that vertex excitations are 
deconfined, and that closed strings without end-points 
have no energy cost. Indeed, taking C = dP gives pre- 
cisely the plaquette operator B p , which commutes with 
the Hamiltonian by construction. A closed vertex type 
string operator that wraps several plaquettes is just the 
product of B p 's for the enclosed plaquettes. However, in- 
terestingly, one can also construct non-contractible ver- 
tex type string operators which flip spins along paths 
C X ,C V which wind completely around the system in the 
x or y direction respectively. The resulting operators 
Wv(C x ),Wv(C v ) flip the parities P y or P x respectively 
and so have the effect of toggling between the four dif- 




FIG. 5. (Color online) This figure shows the different types 
of string operators in the simple 2D lattice models. For the 
toric code, the vertex type string operator simply flips the spin 
on every edge of a path with a x . The plaquette type string 
operator is represented as a line living on the dual lattice, and 
every edge that crosses the line is acted on with a z . For the 
DSem model, the plaquette type string operator is precisely 
the same. The vertex type string operator, however, includes 
additional phases which depend on the spins touching the 
path. A phase (±i)' 1_CT " 2 is associated with the dashed 
(green) edges labelled R which lie just to the right of the 
path; the choice of ± determines the chirality of the string. 
Vertices labelled L are attached to an edge lying on the left 
of the path, and these vertices are associated with a phase 
(— 1) if 1-0 "; )( 1 + CT j ■) w here i/j are the edges on the path just 
before/after the L- vertex. 

ferent ground state sectors in Fig. [4] 

One can similarly define plaquette-type string- 
operators that create pairs of plaquette defects. These 
take the form Wp(C) = Jliec °f' where C is a string 
on the dual lattice as shown in Fig. [5j Vertex defects 
have nontrivial relative statistics with plaquette defects: 
Moving a vertex defect all the way around a plaquette 
defect gives a Berry phase of —1, and vice versa. These 
statistics are a defining element of topological order, and 
will allow us to distinguish the toric code from the dou- 
bled scmion model which we discuss in the next section. 



3. Graphical representation of operators as strings 

Before discussing our next 2D topological lattice 
model, it is useful to review the connection between the 
rules in Fig. [2] and the lattice model. Thus far, we have 
established that the rules in Fig. [2] (b)-(d) can be used 
to evaluate the relative amplitudes of different kets in 
a particular ground state. However, they can also be 
used to depict the action of the string operators defined 
in the previous section. An example is shown in Fig. [6] 
the plaquette operator B p can be represented by a string 
running along the inside of the plaquette p. To deter- 
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FIG. 6. This diagram shows how we can use the fusion rules to 
graphically define a plaquette flipping term (a) using a string 
that lives off lattice shown in (b). The key step occurs in going 
from (b)-(c) where we use Fig.|5|d) to fuse the string into the 
edge. The end result (d) is that all edges on the plaquette are 
flipped, which is precisely the action of Iligap °i ■ 



mine the action of this operator on a particular ket, we 
"fuse" this string with the relevant edge labels of p using 
Fig. id). 

Vertex string operators Jliec a f are defined using a 
similar procedure: we depict the string operator graphi- 
cally as a string running along the path C. The action of 
the operator on a ket is given by fusing this string into 
the edges along C using Fig.[2jd) and (e). (For the Toric 
code the rule (e) is superfluous; however in other models 
it will be necessary as we will assign a definite meaning 
to strings drawn crossing over, as opposed to under, the 
edges in its path (Fig. [7]).) 

Using these graphical representations of states and op- 
erators, it is possible to describe many more 2D topo- 
logical lattice models than the ones we describe here^. 
A careful discussion of the conditions under which these 
graphical procedures give a consistent and unambiguous 
action for the operators in the theory can be found in 
Refs. HandDH 



B. The doubled semion model 

The toric code model shows that it is possible to write 
a simple and exactly solvable spin Hamiltonian exhibit- 
ing topological order. Many generalizations are known 
but here we describe one of the simplest, a second loop 
model known as the doubled semion or DSem model first 
discussed in Refs. [5] and [T3]. The Hamiltonian takes 
the form 

tf = -ELR+E(Il-?) II * (1 ^ )/2 ' ( 4 ) 

v s(v) P i€dp j£s(p) 

B v Bp 




B„W r (C)\<fr) W r (C)B p \^) 




FIG. 7. This diagram illustrates how we can use an off-lattice 
string to define a string-operator W (C). Starting with some 
ket in (a) we define the string operator using the diagram in 
(b), where the black dots denote the end-points of the string 
operator. In order to evaluate this diagram we need to use the 
fusion rule Fig. [2jd) at the point labelled F and the braiding 
rule Fig. [2je) at point X to resolve the string under-crossing 
the lattice edge, (e) Represents a graphical proof that the pla- 
quette operator commutes with the string operator along its 
length; on the left hand side the string acts first, on the right 
hand side the plaquette acts first and we can move between 
these two diagrams by deforming the string locally which is 
allowed as long as string operator does not move through the 
centre of a plaquette. 



where s{p) is the set of six legs radiating from plaquette 
p, as shown in Fig. [Ij and as above s(y) and dp denote 
the three edges entering vertex v, and the six edges bor- 
dering plaquette p respectively. Comparing Eq. Q with 
the toric code Hamiltonian in Eq. ([IJ, we see that B v 
is unchanged but the plaquette terms differ: in addition 
to flipping the edges of p, the new B p operator includes 
a phase which depends on the spins in s(p). Further- 
more, the plaquette part of the DScm Hamiltonian ap- 
pears with a total + sign rather than a — sign. 

For any state \%i oop ) obeying B v \^i oop ) = \^i oop ) (a 
'loop gas state'), it is easy to show that 

[Bp, B v ] \^i oop ) = [B p , B p ,} \f loop ) = (5) 

and furthermore that B p has eigenvalues ±1 on these 
loop gas states. Thus if we restrict our Hilbert space to 
the loop gas states, the Hamiltonian is a sum of com- 
muting projectors and is thus exactly solvable. This fact 
will allow us to solve exactly for the ground states of the 
Hamiltonian, because we will find that these only involve 
closed loop configurations. Though it will not be neces- 
sary for our purposes, with a more complex definition of 
Bp it is possible to write a Hamiltonian that is exactly 
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(a) 



(c) 




lb) 




(d) 



FIG. 8. This figure shows the graphical definition of a B p on 
a ket in the DSem model, (b) shows the definition in terms of 
a string operator. The step (b)-(c) requires the use of a fusion 
rule Fig. [2jd) , so produces a minus sign. Similarly, step (c)- 
(d) requires a fusion, giving a net plus sign. The original ket 
and (d) are related by a single application of B p , therefore 
they appear with opposite sign in the ground state because 
Bp = -1 on the ground state. 



solvable on the full Hilbert space; we direct the interested 
reader to Appendix [Cl 



occur with precisely the same coefficient in the ground 
state (as indicated by Fig. [2jb)). However, the B p = — 1 
condition implies that the two kets involved in 



(-B p ) 




(8) 



Ground States 



occur with a relative minus sign in the ground state su- 
perposition. Therefore spin configurations related by the 
creation or elimination of a closed contractible loop have 
a relative minus sign in the ground state superposition 
(as indicated by the rule in Fig. |2fc)). In the same way 
one can show that configurations related by the fusion 
rule Fig. [2]jd) occur with a relative minus sign. Thus the 
doubled semion graphical rules Fig. |2jb)-(d) give linear 
relations between the coefficients of spin configurations 
in the ground state as illustrated in Fig. [3] 

In the ground state space, the amplitude of any loop 
configuration in can be related, by using the DSem graph- 
ical rules Fig. [2jb)-(d), to the amplitude of one of the 
canonical kets labelled (a)-(d) in Fig. [4] As these new 
rules have minus signs, the ground states are similar to 
those in the toric code except some configurations ap- 
pear with a —1 coefficient, as shown in the superposition 
in Fig. [4| Analogous to the toric code, the ground state 
space is thus split into 2 2 orthogonal sectors distinguished 
by parity operators P x , P y . 



A lower bound on the ground state energy is obtained 
when 



2. Excitations 



B r 



1 for all v 



Br, 



- 1 for all p . 



Any state which satisfies all of these conditions is au- 
tomatically a ground state because it saturates a lower 
bound on the Hamiltonian. Once again such states ex- 
ist. Since the vertex condition is the same, the ground 
states of the doubled semion model are also made of a 
superposition of spin configurations which involve only 
non-intersecting closed loops. However, as we will see 
below, the coefficients in the superposition arc different 
from those of the toric code. 

The condition B p = — 1 implies that the two spin con- 
figurations involved in the equation 



(-B„) 




occur with the same coefficient in the ground state su- 
perposition. Therefore, similar to the toric code, con- 
figurations related by deforming a loop over a plaquette 



(6) The low lying excitations of the doubled semion Hamil- 
tonian appear as pairs of vertex defects where B v = — 1, 
or pairs of plaquette defects where B p = +1. As for the 
Toric code, the operators that create pairs of defects are 
'string operators' acting on all edges along some path C 
joining the pair. Here we will show how to construct 
string operators Wp, (Wv) that violate exactly two pla- 
quettes (vertices) at their end-points, and no other terms 
in the Hamiltonian. It follows that both vertex and pla- 
quette excitations are deconfined: the energy cost of cre- 
ating a pair of either of these is independent of their 
separation. 

The plaquette type string operator is identical to that 
in the toric code, taking the form Wp(C) = FLec a i' 
where C is a string on the dual lattice (see Fig.[5|. Such 
an operator manifestly commutes with every B v , and 
only induces plaquette defects at its endpoints, where it 
anti-commutes with the Y[ ° x part of the B p operators. 
On the other hand, the vertex type string operator from 

(7) the toric code model, Yl ieC erf, fails to commute with the 
modified plaquette operators bordering on C because of 
the v- 1 ~ <T i''' 2 factor in the new definition of B p . To write 
down a vertex type string operator that produces vertex 
violations at its endpoints but commutes with plaquettes 



along its length, Yiiec a f mus t be combined with phases 
depending on the spins on edges touching the path to 
form the following operator Wy (C): 

n>* n i-i)^ 1 -^ 11 ^ n (±*) (i_af)/a - 

«GC k£L vertices /(Ei? legs 

(9) 

To specify the phases such that Wy (C) commutes with 
Bp everywhere, we have had to choose an orientation for 
our string operator (Fig. [5]), and separate the vertices 
into two types: L-vertices are attached to an edge on the 
left of C, while R- vertices are attached to an edge on the 
right of C which we call an R-lcg. Each L-vertex is also 
attached to two edges on C, and these are denoted by 
i, j where i occurs before j on the directed path C. The 
choice of ± determines whether the string has ± 'chiral- 
ity'. The phases in the definitions of the string operators 
imply that two vertex string operators of the same chi- 
rality anti-commute if their paths cross. This indicates 
a sensitivity to the order in which particle pairs are cre- 
ated and destroyed, and can be used to show that vertex 
defects of the same chirality have semionic statistical 
exchanging a pair of these induces a phase of ±i in the 
wave function. 

There is again a simple graphical mnemonic which re- 
produces the required spin flips and configuration depen- 
dent phases along the length of the string away from its 
endpoints. Pairs of 'positive chirality' vertex defects can 
be created by laying a string over the lattice along a pre- 
scribed path and using the rules Fig. [2tb)-(e) to fuse the 
string into the edges; an example of this procedure is 
shown in Fig. [9] The complex conjugate operator cor- 
responds to laying the string under the edges it crosses 
in the lattice, and produces a pair of 'negative chirality' 
vertex defects. 

As in the toric code, it is useful to consider string op- 
erators that do not create excitations (i.e. where the 
curve C has no endpoints) . If C encircles a single plaque- 
tte p, Wy(C) = B p , and hence Wy (C) commutes with 
H. If C traces the perimeter of a cluster of plaquettes, 
Wy(C) = n pG cluster B p an d hence such string opera- 
tors act trivially on the ground states. The operators of 
main interest are again those vertex type string operators 
which follow paths C x , C y which wind completely around 
the system in the x and y direction respectively. The re- 
sulting operators W V (C X ), Wy(C y ) flip the parities P y , 
P x respectively. These non-contractable string operators 
cannot be written as a product of B p operators, but still 
commute with the Hamiltonian and have the effect of 
toggling between the four different ground state sectors 
in Fig. |4j 

We end this section by noting a technical point. The 
above expression Eq. (|9| gives the correct form of vertex 
type string operators away from their endpoints. Captur- 
ing the precise form of the string operator at its endpoints 
requires a careful treatment. The details of this are not 
essential for our purposes, but are shown in Appendix |Cj 



(a) 
W + V (C) 

(c) 




(b) 



JJJ 




(d) 




FIG. 9. This diagram illustrates how we can use an off-lattice 
string to define a positive chirality string-operator W + (C) in 
the doubled semion model; positive chirality corresponds to 
taking +i in the concrete Pauli matrix form of the operator 
in Eq. Q. Starting with some ket in (a) we define the string 
operator using the diagram in (b), where the black dots denote 
the end-points of the string operator. In order to evaluate this 
diagram we need to use the braiding rule Fig.[5Je) at the string 
over-crossing labelled X, yielding a factor of i. We then use 
and the fusion rule Fig. [2jd) at point F to fuse the string into 
the edge, yielding an additional factor of —1. 



III. SIMPLE 3D TOPOLOGICAL LATTICE 
MODELS ON CLOSED MANIFOLDS 

Levin and WerP suggested a method generalizing a lim- 
ited subset of their exactly solvable models to 3D. Their 
method can be used to generalize the 2D toric code to 
3D but cannot be used to generalize the DSem model to 
3D. In this section we follow the different prescription 
of Walker and Wang^ which allows us to generalize both 
the 2D toric code and DSem models to 3D. The resulting 
3D lattice models are described by a two-state system 
o~ z = zbl on each edge of the trivalent lattice shown in 



Fig. 10 T he gen eralization of the 2D toric code to 3D is 
well knowrpSEfil. We nevertheless begin by reviewing the 
key signatures of topological order in this model (ground 
state degeneracy, possible defects and their statistics) for 
later contrast with the properties of the 3D cousin of the 
DSem model. In Sec. |III A| we discuss the ground state 
degeneracy of the 3D Toric code on the 3-torus, and de- 
scribe the low-lying excitations which take the form of de- 
confined point defects and linearly confined vortex rings. 

In Sec. |III B| we present the 3D analogue of the 2D 
DSem model (the '3D semion model'), which is our first 
example of a confined WW model. Again we will charac- 
terize the phase by studying its ground states and excita- 
tions. We find that the ground state is non-degenerate on 



the 3-torus, and in Sec. IIIB 3 we relate this to the fact 
that its point defects are confined in the bulk. Combined 
with the study of these models on a manifold with bound- 
ary which we present in Sec. |IV[ this section provides a 
firm foundation for understanding the effective field the- 
ories for the lattice models in Sec. [V] Furthermore, the 
3D toric code and semion models capture much of the 
qualitative behavior of more general Walker- Wang mod- 
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FIG. 10. This figure shows the point splitting and fixed pro- 
jection used to define the 3D lattice for the Walker-Wang 
models. The dot in the middle of each bond represents a spin 
variable. 



els which we study in Sec. VI 



A. The 3D toric code 

V iewed as a Walker-Wang model, the 3D toric 
cod e! 13 * 16 -* Hilbcrt space consists of a two state system 
er z = ±1 on each edge of the lattice shown in Fig. [10] 
The Hamiltonian takes the form 



H 



z n <-£n 



(10) 



where s(v) is the set of three edges attached to vertex v 
and dp is the set of ten edges of a plaquette p (bold edges 
shown on the 3 types of plaquette in Fig. [l4|a)-(c)). As 
in 2D, the B v take the values ±1 depending on whether 
there are an even/odd number of down spins on the edges 
coming into vertex v, and B p flips the spins on each edge 
of p. The fact that B p flips a pair of spins at vertex v 
implies that \B P ,B V \ = 0, and once again [B V ,B V >] = 
[Bp, B p i] = 0, so the model is exactly solvable. 



Ground states of toric code on T 



The ground state space is defined by the conditions 
Bp = B v = 1 for all vertices and plaquettes. As in the 
2D case in Sec. |IIB T] the condition B v = 1 forces the 3D 
ground state to be a superposition of closed loops. Re- 
call that in 2D, the condition B p = +1 implied that the 
rules in Fig. [2] (b)-(d) relate the amplitudes of different 
spin configurations. Exactly the same type of calculation 
shows the rules in Fig. [2] (b)-(d) relate the amplitudes of 



spin configurations in the ground state space of the 3D 
toric code. 

In the 2D toric code on the 2-torus, the amplitude 
of any spin configuration in the ground state could be 
related to the amplitude of one of four canonical con- 
figurations (shown in Fig. [4]) resulting in four degener- 
ate ground states. Analogously, on the 3-torus one can 
show that any configuration of closed loops appears in 
the ground state with the same amplitude as one of the 
eight canonical configurations shown in Fig. |11| These 
2 3 configurations can be labelled by three parities 



p nx = n 



x,y,z 



(11) 



which take values ±1 depending on whether an even or 
odd number of loops wind around the n-cycle of the torus. 
Here z± is the set of all edges emanating from the plane 
(z = 0) in the -(-indirection, and similarly for x±,y±. If 
we take an equal superposition of all configurations re- 
lated to one of the eight canonical configurations by the 
equivalences in Fig. [2ft))- (d) (as shown in Fig.[l2| then it 
is easy to verify that we get a ground state. This shows 
us that the ground state degeneracy is 2 3 , and each dis- 
tinct ground state is labelled by the three eigenvalues 
P x± , P y± , P z± = ±1. On more general closed manifolds 
the ground state degeneracy is 2 bl , where b\ is the num- 
ber of independent non-contractible cycles on the mani- 
fold (also known as the first Bctti number). 



2. Excitations in the toric code 



The Hamiltonian ( 10 ) has two types of excitations: 
pairs of vertex defects where B v — — 1, and lines of pla- 
quette defects where B p = —1. We can create a pair of 
vertex defects with an operator 



W V (C 



AB 



n * 



(12) 



where Cab is a path connecting the positions of the de- 
fects A and B. Graphically we represent the string oper- 
ator by laying a string along Cab , where it is understood 
that the operator acts on kets by fusing this string into 
the edges using the rules Fig. [2ft))- (e). The operator 
commutes with the Hamiltonian except at its endpoints, 
and so the defects are deconfined. 

As in 2D, a closed vertex type string operator that 
wraps around the boundary of several plaquettes is 
just the product of B p s for the enclosed plaquettes 
and so trivially commutes with the Hamiltonian. Non- 
contractible string operators (i.e., string operators that 
wrap around the periodic boundary conditions) are more 
interesting because they commute with the Hamiltonian 
but cannot be expressed as a product of Bp's. Oper- 
ators of this form toggle between the different ground 
state sectors discussed in Sec. |III AT] for example, a non- 
contractible string operator wrapping the z-direction of 
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FIG. 11. Shown are the eight canonical configurations with 
periodic boundary conditions in all three directions, where the 
thick black lines indicate edges with <r z = — 1. The underlying 
lattice is not drawn, for simplicity. Any basis ket without 
vertex violations can be related to one of these eight by using 
the graphical rules in Fig. [2] In the case of the toric code, the 
ground state splits into eight orthogonal sectors labelled by 
these eight kets; two of these states are shown in Fig. [12] 



the torus toggles between the ground states in Fig. 12 'a) 
and (b). 

Thus, vertex defects in 3D are much the same as they 
were in 2D, being thought of as the end-points of a string 
operator. Plaquette defects, on the other hand, behave 
quite differently in 2 and 3 dimensions. They no longer 
appear at the end-points of string operators, but rather at 
the boundary of surface operators. To create a plaquette 
defect, one acts on an edge with a a z . This, however, 



creates defects in all four plaquettes associated with the 
edge. More generally, if we pick a surface on the dual 
lattice and act with a z on each edge cutting the surface, 
then the resulting operator 



(13) 



ies 



creates plaquette defects along the boundary dS of the 
surface S (see Fig. [l3|a)). This is because a z acts on an 
even number of edges of each plaquette cutting S, except 
for those lying on the boundary dS which only have a z 
acting on one of their edges. Thus this type of surface 
operator has an energy cost which scales linearly with 
the length of the boundary dS of the surface. 

One of the defining features of the topological order of 
the 3D Toric code is the mutual statistics between the 
point-like vertex defects and the vortex lines. Moving 
a vertex defect all the way around a line defect gives 
a Berry phase of —1. The exchange process is defined 
as follows. First act on a ground state with the two 
operators to form Wy (Cab)Wp(S) | GS) as shown in 
Fig. 



13 



a.). We now take the point defect B and thread 
it through dS, and then annihilate A with B. This results 
in a string operator Wv{C c i OS ed) enclosing the plaquette 
defects on dS as shown in Fig.[l3^b). The string operator 
encircling the line defect can be written as a product of 
plaquette operators Wv(C c i OS ed) = Yipen Bp, where TZ is 
any surface bounded by C c i ose d- However, one of the B p 
operators in TZ will lie on the line defect, and therefore 
take value —1. Thus a full exchange between the two 
types of defects leads to a sign of —1. 



B. The 3D semion model 

We now define the '3D semion model', our first exam- 
ple of a Walker- Wang model which has not previously 
been examined in the literature. To define the operators 
in this model, we first fix a projection of the trivalent 

(Note that we will 
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lattice onto 2D, as shown in Fig. 
always assume the lattice is defined on an orientable man- 
ifold.) The Hamiltonian takes the form 



H = - 




(14) 



p i€dp j£s(p) 



where n = ~ (1 — a z ) and, as for the 3D toric code, s(v) 
is the set of three legs attached to vertex v, dp is the set 
of ten edges bounding plaquette p, and it includes the 
two privileged blue and two privileged red edges used in 
the definition of B p above, while s(p) is the set of ten 
edges radiating from plaquette p (see Fig. |14[a)-(c)). 
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FIG. 12. This figure shows two examples of ground states of the 3D toric code, (a) consists of an equal superposition of all 
kets related to Fig. |llfr , by Fig. [5]^a)-(d), and has P x± = P y± = P z± — 1. (b) consists of an equal superposition of all kets 
related to Fig.JTTJd) by the same graphical rules, and has P x± = P y± = 1 but P z± = — 1. There are a further six ground states 
corresponding to the other kets in Fig. |11| The underlying lattice is not drawn, for simplicity. 



As in 2D, for states |*j oop ) obeying B v \^ loop ) = |*; oop ) 
(so-called 'loop gas states'), we can show that 

[B p , B v ] \f loop ) = [B P1 By] |* ioop ) = , (15) 

and furthermore that B p has eigenvalues ±1. Thus, if 
we restrict our Hilbert space to the loop gas states, the 
Hamiltonian is a sum of commuting operators. Once 
again, this fact will allow us to solve exactly for the 
ground state of the Hamiltonian, because we will find 
that the ground state is a loop gas state. We can 
also define B p by using the graphical rules: Thread a 
string around the boundary of plaquette p as shown in 
Fig.[l4|d)-(f), and fuse it into the edges using the graphi- 
cal rules (Fig.[2jd) and (e)). The string will have a choice 
of under-crossing or over-crossing two edges in s(p); these 
are the U and O edges marked in Fig. |l4[d)-(e). In 
the Walker- Wang prescription the plaquette string will 
under-cross the U edge, and over-cross the O edge, a 
choice which ensures adjacent plaquettes commute. Note 
that in 3D we need to use the crossing rule Fig. J2je) in 
order to define the plaquette operators graphically, while 
in 2D we only needed this rule to define the string oper- 
ators. 



Unlike the Toric code Hamiltonian (10), the complex 



phase in ( 14 ) ensure that this Hamiltonian explicitly 
breaks time-reversal symmetry. However, it is natural 
to define the parity operator such that it also acts by 
complex conjugation on H (which corresponds to choos- 
ing an opposite-parity projection of the lattice). With 
this definition PT is a symmetry of the semion model. 



1. Ground state of 3D semion model on T 

We first contrast the topological order in 3DSem with 
that of the 3D toric code by comparing the ground states 
of the two models. The 3D toric code has a ground state 



degeneracy of 2 3 on the 3-torus, where each ground state 
is labelled by three parity eigenvalues P n± — ±1 with 
n = x,y,z. We will find that 3DSem has precisely one 
ground state on the 3-torus, with all three parity values 
P n± = +1. 

In analogy with the 2D DSem model, the 3D semion 
ground state is defined by the conditions B v — 1 and 
B p = — 1 for all vertices and plaquettes. The vertex 
condition implies that the ground state is a loop gas. In 
the 2D DSem model, the condition B p = — 1 implied 
that the rules in Fig. [2] (b)-(d) related the coefficients of 
spin configurations. In the 3D semion model, something 
similar happens: the rules Fig. [2] (b)-(d), as well as (e) 
relate the amplitudes of different configurations in the 
ground state superposition. To see the rule Fig. ge) m 
action, note the equation 




, (16) 



along with the condition B p = — 1, implies that the 
ground state amplitudes of these two configurations must 
be related^ by a factor of — i. This is in accordance with 
the graphical rules, which tell us that 




P 



(17) 



In other words, a twist in a loop segment (like that on 
the left hand side of Eq. (16l) should be associated with 




FIG. 13. (Color online) This diagram illustrates the exchange 
between a line defect created by Wp(S) and a pair of point 
defects created by Wv{Cab), both shown (a) for the 3D toric 
code. We show the process on a cubic lattice for simplic- 
ity, but the same reasoning applies on the trivalent Walker- 
Wang lattice. The defect B is threaded through the surface 
5 bounded by the line defect, and annihilated with A. The 
resulting closed loop can be written as a product of plaquette 
defects, yielding a result of —1 because one of the plaquettes 
lies on the line defect. 
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FIG. 14. (Color online) (a)-(c) show the three different types 
of plaquette occurring on the lattice. The edges have been 
made either black, bold black, bold blue (labelled B) or bold 
red (labelled R) to aid in the definition of the 3D semion 
plaquette operator in Eq. Q. The set dp contains the ten 
bold edges, while the set s(p) contains the ten black edges. 
(d)-(f) show how to define B p operators in terms of a string 
picture, where it is understood that the string is fused into 
the edges using Fig.[5|d) and (e); notice that the string under- 
crosses the edge labelled U, but over-crosses the edge labelled 
O. 



the plane perpendicular to the n-direction. This follows 
immediately from the identity 



a phase of —i relative to a straight loop segment (such 
as the one on the right of Eq. (16)). With the aid of 
these graphical rules, we show that the 3D semion model 
has a non-degenerate ground state on the 3-torus. In 
the case of the toric code, there were precisely 2 3 ground 
states; each ground state was labelled by the three par- 
ity eigenvalues {P n± , where n = x, y or z), and was 
formed by making an equal amplitude superposition of 
all configurations related to a canonical ket (with the cor- 
responding parity) by the graphical rules. In the case of 
the 3D semion model, we can again use local rules to 
relate any loop configuration to one of the eight canon- 
ical configurations in Fig. |11[ but this only proves that 
there is a ground state degeneracy of at most eight. In 
fact, there is only a single ground state, which has all 
P n± = 1; it is a superposition of all loop configurations 
related to Fig. fTTVa) , with relative phases given by the 
rules Fig. [2]jb)-(e) as shown in Fig. 15 



What happens to the other configurations, which have 
some P n± = — 1? None of these configurations occur in a 
ground state superposition because they are all excited: 
P n , = —1 implies the existence of plaquette defects in 



P, 



n 

pen„ 



(18) 



which we prove graphically in Fig. |16[ where H n± is the 
set of plaquettes lying in a plane perpendicular to the n 
direction. We see that any negative parity P n± = — 1 is 
incompatible with the ground state conditions that B p = 

— 1 for all p. Technically we still need to prove that the 
ket Fig. 11 'a) has an overlap with the ground state. This 
can be shown by explicitly expanding the ground state 
projector V, as we do in Appendix [FJ 

Although we have worked specifically with the 3- 
torus, the above method appears to generalize to any 
(orientable) manifold without boundary. Hence, as 
promised, we have shown that the ground state of 3DSem 
is unique on any (orientable) manifold without boundary 

- unlike both the 3D Toric code, and the 2D semion 
model. 
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FIG. 15. This figure shows the unique ground state of the 3D semion model on the 3-torus. It is a superposition of all 
configurations with all three parities P x± — P y± — P z± — 1. To calculate the coefficient of any ket use Fig. [2] For example, a 
ket with a single loop has coefficient —1 using Fig.[2jc). One can show that the ket with the twisted loop has coefficient —i by 
using Fig. [2|e) then (c). 



2. Bulk excitations in the 3D semion model 

A second important characterization of topological or- 
der is the identification of the low lying excitations and 
their statistics. In the 3D toric code, the low lying exci- 
tations are vertex string operators which produce decon- 
fined point defects, and sheet operators which produce 
a line of plaquette defects at their boundary. The 3D 
semion model is very different because its point defects 
are confined: We will show that any attempt to separate 
a pair of point defects leads to the creation of a line of 
plaquette violations B p — 1 connecting them. 

First we note that, unlike in the 3D toric code, flip- 
ping spins along a path using Ili^CAe °^ does not pro- 
duce deconfined point excitations. This operator fails 
to commute with the B p along its length due to the 
(T z -dependent phase terms in Eq. (14). The same kind 
of issue arose in 2D in the DSem model (Sec. II B 2 1, 
where we found that the string operators producing de- 
confined point defects were not simply strings of a x op- 
erators; in addition to flipping the spins along a path, 
the DSem string operators included compensating phases 
(depending on the states of edges touching the path) 
which ensured they commuted with the plaquettes along 
their length. However, in 3D there is no such assign- 
ment of phases, a fact we will prove in the next section 
Sec. |IIIB 3| For now, we will demonstrate how a particu- 
larly important type of string operator fails to commute 
with plaquettes along its length. 

In the 2D DSem model there was a neat graphical rep- 
resentation for the action of string operators (Sec. II B 2 1: 



lay a string along a path either over or under the lattice, 
and use the graphical rules Fig. [2]jb)-(e) to fuse it into 
the edges. Now attempt to use this prescription in 3D: 
lay a string along a path C connecting two defects, taking 
note of which edges the path over and under-crosses, and 
fuse it into edges using the rules Fig. [2]jb)-(e) to form 
an operator Wy(C). In addition to flipping spins along 
the path, the resulting operator includes phases which 
ensure it commutes with many of the plaquettes touch- 
ing C. How ever, if the path C threads a plaquette p, as 
in Fig. 



17 



then Wv(C) flips the eigenvalue from the 
low energy state B p = — 1 to an excited state B p = +1; 
we prove this statement in the caption to Fig. |17| Thus, 
attempting to create a pair of point defects with Wy[C) 



leads to a line of plaquette defects. 

We might have expected from the outset that point 
particles are confined in the 3D semion model. In the 2D 
DSem model, point defects had relative semionic statis- 
tics: they were neither fermions nor bosons. The 3D 
semion model is built on the same semionic graphical 
rules, so we might have expected point defects in 3D to 
behave like semions. However, in 3D, free point-like par- 
ticles are either fermions or bosons. The model complies 
with this requirement by confining the pairs of semionic 
point particles. Intriguingly, we will later find in Sec. |IV| 
that deconfined semionic defects do exist in these 3D lat- 
tice models, with the caveat that they are pinned to a 
boundary. 



3. Connection between confinement and non- degeneracy 

We saw above that the toric code has degeneracy 2 3 
on the 3-torus. This degeneracy is associated with the 
existence of de-confined point excitations in the follow- 
ing way. If the point defects are de-confined this means 
they are associated with the endpoints of vertex type 
string operators which commute with the Hamiltonian, 
except at their endpoints. Therefore closed string op- 
erators that wrap around the non-contractible cycles of 
the 3-torus have no energy cost. It is easy see how 
these operators toggle between the different sectors of the 
ground state. For example, acting on the ground state in 
Fig. fT2|a) with the operator Iljec- a j yields the ground 
stateFig. [l2^b), where C z is any path that goes once 
around the z-cycle of the torus, and none of the other 
cycles. This should be contrasted with the 3D semion 
model, which has a single ground state on the 3-torus. 
In this unique ground state (Fig. [l5^a)), all parities are 
trivial (P n± — !)• This observation proves our previous 
suggestion that point defects are confined in the bulk 
of the 3D semion model: if they were not confined, we 
would be able to form a closed string operator with no 
energy cost, and use it to thread a single string around a 
non-contractible cycle in the torus, forming another dis- 
tinct ground state with non-trivial parity. This would 
contradict the fact that there is only one ground state. 

Unlike the 3D Toric code, then, the 3D semion model 
has no topological order in the conventional sense. Com- 
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FIG. 16. This figure shows the graphical procedure for eval- 
uating the operator Ilpgn B p on ^ ne 3-torus. \l z± denotes 
the set of plaquettes in the plane perpendicular to the z di- 
rection, while z± represents a set of edges encircled by loops 
in (c); \H Z± \ and \z±\ represents the sizes of these sets. To 
begin, we evaluate (b) by fusing each neighboring string. This 
yields the picture in (c), where loops are left encircling edges 
coming up out of the plane as well as a factor of (— 1)^ 
for all the fusions. The loops encircling the vertical legs 
have the same action as — a z , as shown in (d). This gives 
n p6 n z , B p = (-l) / ~' Zil Y\ iGz o-j. Careful counting shows 

(-1)'- 



= (-1)1 



I , which gives Eq. 
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mensurate with this, its bulk entanglement entropy also 
vanishes, as we show in Appendix |D| 

Having established that the 3D semion model is topo- 
logically trivial on manifolds without boundary, we now 
turn to study the model on a manifold with boundary, 
where we will find that its physics is markedly different. 



FIG. 17. (Color online) This diagram shows how a vertex type 
string operator Wv(C) in the 3D semion model violates any 
plaquette p it threads. Using the graphical rule in Fig. [2^e) 
shows that the state B p Wv{C) | GS) in (a) is equal to the 
state -W V (C)B P \ GS) = W V (C) \ GS) in (b). Hence B p = 
+ 1 on the state Wvifi) \ GS), and so the string operator has 
excited the plaquette p out of its low energy state B p — — 1 
to an excited state B„ = +1. 



IV. 3D TOPOLOGICAL LATTICE MODELS ON 
MANIFOLDS WITH BOUNDARY 



In the previous section we investigated the 3D semion 
model on manifolds without boundary; all excitations 
were confined and we connected this with the unique- 
ness of the ground state. In this section we investigate 
the 3D semion model on manifolds with boundary, find- 
ing that the presence of the boundary leads to important 
qualitative changes in the topological order: There exist 
deconfined chiral anyonic excitations on the surface of the 
manifold, and these are associated with the existence of 
multiple ground states. In contrast, the presence of the 
surface will have little qualitative effect on the topologi- 
cal order in the 3D toric code; in particular, there are no 
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FIG. 18. (a) Shows a lattice representing a solid donut D 2 x S 1 
where periodic boundary conditions are enforced in the z- 
direction. Following that, (b)-(d) show how to define various 
string operators by laying a string over the lattice and fusing 
it into the edges using graphical rules, (b) Shows the action of 
a closed non-contractible string operator W(C Z ) that toggles 
between the two ground states, which are distinguished by 
whether an odd or even amount of flux goes around the z- 
cycle. (c) Defines the operator W(Cj,), which also commutes 
with the Hamiltonian but anti-commutes with W(C Z ). (d) 
Shows how to define deconfined surface excitations at A and 
B using an open string operator W(Cab)- 



new excitations. This section is organized as follows. In 
Sec. |IV A| we define the boundary conditions for the toric 
code on the solid donut, and then describe its ground 
states and excitations. Then in Sec. IIVBI we define the 
analogous boundary conditions for the 3D semion model, 
and show that there are two ground states which we tie 
to the novel deconfined anyonic surface excitations. 

The choices of boundary conditions in 3D are similar 
to those in 2D, which have been discussed at length in 
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FIG. 19. This figure shows an example of the surface of a lat- 
tice that has been 'smoothly' cut-off. The boundary Hamil- 
tonians of both toric code and 3D semion models will include 
vertex and plaquette projectors for every (trivalent) vertex 
and plaquette in the surface. 



Refs. [18] and [T9j In this section we focus on one of the 
3D analogues of this, the 'smooth' boundary shown in 
(We will also briefly discuss the rough bound- 



ary condition, which does not admit an exactly solvable 
surface Hamiltonian, in Appendix [E]) We will show that 
adding a chiral 2D anyon layer to this boundary can de- 
stroy the topological order associated with the surface. 



The toric code on the solid donut 



In this section we examine the toric code on the solid 
donut manifold D 2 x S 1 , where D 2 is a 2D disk and S 1 
is a circle. The choices of boundary conditions in 3D 
are similar to those in 2D, which have been discussed at 
length in Refs. [TBI and IT5] . Here we will consider one of 
the 3D analogues of this, the 'smooth' boundary shown 
for a small (one plaquette thick) version of the solid donut 
in Fig. [18J and more generally in Fig. [T9| The toric code 
Hamiltonian is again defined by Eq. (10), where v runs 



over all trivalent vertices, and p runs over all plaquettes, 
including those plaquettes and vertices lying on the sur- 
face of the manifold. These boundary conditions gap the 
surface to vertex and plaquette defects, which will ensure 
that the ground state degeneracy is independent of the 
number of sites on the lattice. 

The ground state is defined by B v = 1 and B p = 1. 
This again implies that the ground state is a loop gas 
state, and that we can relate different configurations us- 
ing the graphical rules Fig. [2|a)-(d) as before. In the 
case of the solid donut, there are two 'canonical' config- 
urations to which all other loop gas states can be de- 
formed, each distinguished by the winding parity P z± . 
There are 2 ground states labelled by P z± = ±1, and each 
is an equal sum superposition of all loop configurations 
related by graphical rules to the corresponding canoni- 
cal configuration. The ground state degeneracy is 2 bl on 
more general manifolds with these boundary conditions, 
where &i is the number of independent non-contractible 
cycles. 
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1 . Excitations 

The excitations of the toric code on the solid donut 
are much the same as they were for the 3-torus. There 
are deconfmed point particles in the bulk and on the sur- 
face, created by vertex type string operators Wy{C) — 
Yljec a j- ^ s before, we are particularly interested in 
string operators which are closed and non-contractible 
like Wv{C z ), where C z is any path that wraps once 
around the (non-contractible) z-cycle of the solid donut. 
This operator commutes with the Hamiltonian, but flips 
P z± , so toggles between the two possible ground states. 

As on the 3-torus, we can also use operators of the form 
Wp(S) = n. ieS &i to create lines of plaquette defects 
along the boundary dS of a surface on the dual lattice S. 
An interesting feature of our chosen boundary conditions 
is that Wp(S) does not have an energy cost for parts of 
dS which lie on the surface. We can use this to create 
open lines of plaquette defects in the bulk, provided they 
end on plaquette defects on the surface of the manifold. 
To create a line of plaquette defects connecting surface 
plaquettes p\ , pi along a line C in the bulk, find a surface 
5 which has C as part of its boundary. Make the rest 
of the boundary of S lie on the surface of the manifold. 
The resulting operator Wp(S) produces the line C of 
plaquette defects ending at pi,P2, but no other defects. 

In summary, the presence of the surface does not 
change the nature of defects in the 3D Toric code; they 
still consist of deconfmed vertex defects, and linearly con- 
fined vortex lines (although these lines can now be open 
provided they end on the boundary). The excitations 
are either deconfmed vertex defect, or lines of plaquette 
defects. This is very different to the behavior of 3DSem 
model on the solid donut, which we study next. 



B. The 3D semion on the solid donut 

Evidently, the presence of a 'smooth' surface has lit- 
tle effect on the toric code; the ground state degeneracy 
still only depends on the number of independent non- 
contractible cycles, and we can toggle between all sectors 
of the ground state by using vertex string operators that 
lie in the bulk or on the surface. We now contrast this 
with the 3D semion model, using the same 'smooth' cut- 
off of the lattice and analogous boundary conditions. On 
the solid donut, the 3D semion model has Hamiltonian 



B„ in Eq. ( 14 ) . It is easy to show that the new plaquette 



* = E 



B„ 



E 



nsuri 



(19) 



p in bulk 



p on surface 



where v runs over all trivalent vertices, including those 
on the surface. For those plaquettes p with both dp and 
s(p) in the bulk, the definition of plaquette operator B p 
is as it was in Eq. (14). A plaquette p lying near or 



entirely in the surface will have some edges missing from 
s(p); the corresponding operator £?p Urf is defined simply 
by omitting these edges from the original definition of 



operators still satisfy (B p url ) 2 — 1, and therefore have 
eigenvalues ±1. These boundary conditions are again 
chosen to gap the surface to vertex and plaquette defects. 

Having defined the operators of the Hamiltonian, we 
note that the ground state is defined by B p — —1 and 
B v = 1. These conditions again imply that the ground 
state is a loop gas state, and we can use the graphical 
rules of Fig. [2ja)-(e) to relate the coefficients of different 
configurations. In the present case there are two 'canon- 
ical' states (labelled by P z± = ±1) to which all other 
loop configurations can be related using the graphical 
rules. This suggests that the ground state degeneracy 
is at most two. In fact it is precisely two, and one can 
readily prove this using the methods of Appendix [F] 



1. Surface excitations in the 3D semion model on a 
manifold with boundary 

So far we have shown (See also Appendix [f]) that the 
toric code and 3DSem have degeneracy 2 on the solid 
donut. However, they exhibit very different kinds of 
topological order, a point we clarify by comparing their 
excitations. In the 3D toric code, point defects are de- 
confined on the surface and in the bulk. Being deconfmed 
3D particles, they must be bosonic or fermionic; it is easy 
to show that they are bosonic. In this section, we show 
that 3DSem has deconfmed vertex defects on the sur- 
face of the manifold. These defects are chiral semions. 
However, we cannot have semions in 3D, only bosons or 
fermions! The model complies with this requirement be- 
cause the bulk confinement of vertex defects (Sec. |IIIB~2 ) 
pins the deconfmed defects to the surface. Thus they are 
effectively 2D particles, and therefore are allowed to have 
anyonic statistics. 

In 2D we found that there were two chiralities of de- 
confined point defects; positive chirality and negative chi- 
rality corresponded to the choice of +i or —i in Eq. (|J). 
The string operators creating the point defects had a neat 
graphical mnemonic demonstrated in Fig. [7] To act with 
a string operator, lay a string along the prescribed path 
and fuse it into the edges using the graphical rules Fig. [2] 
Positive chirality strings lie above the lattice, while neg- 
ative ones lie below. 

However, on the boundary of 3-manifolds there is only 
one species of deconfmed vertex defect. To create a pair 
of these vertex defects at A and B, lay a string just out- 
side the surface as shown in Fig. 18 b) along a path Cab 



connecting A to B, and fuse it into the edges using the 
graphical rules. (Recall that laying the string outside the 
surface, as opposed to inside, corresponds to a particu- 
lar choice of phases ±i when this string crosses a spin- 
down edge). Such a string operator (W(Cab)) has phase 
choices which ensure it commutes with the plaquettes 
along its length; graphically, this corresponds to the fact 
that the strings do not thread any plaquettes. In con- 
trast, strings lying just beneath the surface necessarily 
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produce confined point excitations because they excite 
the plaquettes along their length in the same way that 
bulk string operators do (see Fig. 17). This preference 
for one chirality over the other is a symptom of the fact 
that the model breaks time reversal symmetry because 
of the i phases appearing in the Hamiltonian Eq. (14). 



The chiral phase choice (needed in defining the surface 
string operators) is such that two vertex string operators 
on the surface will anti-commute if their paths cross (see 
Sec. Ill B 2 ). This indicates a sensitivity to the order in 



which particle pairs are created and destroyed on the sur- 
face, and can be used to show that the deconfined surface 
vertex defects have semionic statistics^. 

These chiral surface vertex defects are deconfined be- 
cause the surface strings used to create a pair of defects 
commute with the plaquettes along their length. There- 
fore if we form a non-contractible surface string Wv(C z ) 
which winds once around the z-cycle of the torus as 



shown in Fig. 18 b), it will commute with the Hamilto- 
nian and flip the sign of the parity operator P z± . This op- 
erator toggles between the two ground states mentioned 
in the previous section. 

For the 3D semion model on the 3-torus, we found 
there was a single ground state and that any state with 
(for example) P Zl = —1 is excited. This latter fact fol- 
lowed from Eq. (18), which showed that P z± — — 1 im- 



plies the existence of plaquette defects in the z± plane. 
What is the analogue of Eq. ( 18 1 on the solid donut, and 
how does it allow a ground state with P z = — 1? Using 



the same methods used to derive Eq. (18), we can show 
that 



P z± = -W(Ct) ]"[ (-B p ) 



(20) 



p G 2j_ plane 



where WiC^) is another string operator on the sur- 
face which commutes with the Hamiltonian, defined in 
Fig. [l8jc). We see immediately that P z± = — 1 is not 
incompatible with the ground state condition B p = — 1, 
provided W(C<p) — 1. Note that, when there are no pla- 
quette defects, P z± = —W(C$), and so Wv{C^,) anti- 
commutes with Wv{C z ). 



2. Other boundary conditions 



the surface topological order is to our choice of boundary 
conditions? One might imagine that the surface topolog- 
ical order might be altered or even destroyed by adding 
some 2D model to the boundary, i.e., by choosing an ap- 
propriate perturbation to the surface of the 3D system. 

To completely destroy the surface topological order, let 
us begin with a smooth boundary, and add a Laughlin 
v = 1/2 quantum Hall system of opposite chirality to 
the surface. We take the quantum Hall system to have 
a finite density of quasi-particles, which are not strongly 
pinned in space. Next, we perturb the surface Hamil- 
tonian by adding a very strong attraction between the 
vertex defects on the surface of the 3D model and the 
quasi-particles of the quantum Hall system in such a way 
that each quasiparticle will bind exactly one surface ver- 
tex defect (we must make the binding sufficiently strong 
so as to overcome the energy gap to creating vertex de- 
fects, but much less strong than the quantum Hall gap). 
The key to altering the surface topological order is that 
these composite particles have bosonic statistics with the 
semionic braiding phases of the vertex defects precisely 
canceling the phases from the quantum Hall quasiparticle 
braiding. 

We now introduce a small kinetic term to the Hamil- 
tonian so as to allow the vertex defects to hop from site 
to site. (Such a term can indeed be written down so as 
to commute with all of the bulk plaquettes. Details of 
this are given in Appendix D). We therefore expect that 
a bose condensate of these excitations will form at the 
surface. Interestingly, this bose condensate confines the 
chiral anyons present in the smooth boundary, as these 
braid non-trivially with the condensate^. This confine- 
ment ensures that the ground-state degeneracy is lifted, 
leaving only the even parity ground state (P z± — 1). 

In summary, unlike in the 3D toric code, the 3DSem 
model behaves very differently on systems with and with- 
out boundaries. In closed systems it appears topologi- 
cally trivial (unique ground state, confined excitations). 
With the an exactly solvable "fixed point" boundary, 
however, we find deconfined chiral semions living at the 
boundary, and a corresponding ground-state degeneracy 
when the boundary is topologically non-trivial. This 
topological order, which arises purely at the boundary, 
can nonetheless be destroyed by adding a sufficiently 
strong perturbation to the surface. 



We have seen that in the presence of a smooth bound- 
ary, the confined Walker- Wang models are topologically 
ordered, with a topological ground state degeneracy and 
deconfined chiral anyons on their boundaries. An exactly 
solvable zero correlation length (fixed point )P model ca- 
pable of describing chiral topological order is striking: 
as discussed in the introduction, there is good reason to 
believe that such a description is only possible at the 
boundary of a 3D system. 

However, it is certainly possible to generate chiral 
anyon theories in 2D lattice models that are not fixed 
point Hamiltonians. Given this, we might ask how robust 



V. EFFECTIVE FIELD THEORIES FOR THE 
3D TORIC CODE AND 3D SEMION MODEL 

Having described the ground states of the 3D toric code 
and 3D semion models, as well as the properties of their 
excitations, we are ready to investigate effective field the- 
ories for the two lattice models. In Sec. IV Al we review 
(3 + 1)D abelian U(l) bF theory and its connections to 
the Toric code^, and show that it captures two impor- 
tant aspects of topological order in the 3D toric code: 
the ground state degeneracy, as well as the types of ex- 
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citations and their statistics. Then in Sec. VB we show 



that adding a ,pa 'b p V b p(J term to the bF action simi- 
larly captures important aspects of the topological order 
of the 3D semion model. In particular, as outlined in 
[5], this field theory gives rise to a chiral Chern-Simons 
anyon theory on the surface of the manifold, which ac- 
curately describes the deconfined semionic vertex defects 
found on the surface of the lattice model in Sec. IIVB II 
We then go further to show that the field theory forbids 
pairs of free bulk point particles by insisting that they are 
bound to to the ends of line-like objects; this corresponds 
to the fact that vertex defects in the bulk of the confined 
Walker- Wang lattice model arc found at the ends of lines 
of plaquettes defects (Sec. IIIB2). 



A. The bF description of the toric code 

In this section we review the well-known fact that 
(3+l)D abelian bF theorjPH captures the topological be- 
havior of the 3D toric code^. bF theory has two fields: A 
vector field Ap , and an antisymmetric rank-2 tensor field 
bp V . In what follows we will work, for sake of simplicity, 
in 4D Euclidean space with cartesian co-ordinates, where 
there is no distinction between upper and lower indices 
because the metric takes the form gp V — S^v The field 
theory has an action 

5 bF [J,S] = Jd 4 x(J^e^bp, v d p A a + ^^+A^jA , 

(21) 

where J p and E A " / source the world-lines of point parti- 
cles, and the world-sheets of line-like objects respectively. 
More precisely, J p is defined as 



J"(a;) = / daS {4 \x- A») 



dX p 
da 



(22) 



where {X p (a)} is a particle world-line parameterized by 
a. Similarly, is defined as 



On a more general manifold 23 A4, it is 2 bli - M \ where 
b\ (M) is the number of independent non-contractible 
cycles in the spatial manifold. This matches our calcu- 
lation of the lattice model ground state degeneracy in 
Appendix |F 2| In addition to giving the correct ground 
state degeneracy, the (3+l)D bF theory also yields the 
correct statistics between point particles (sourced by J p ) 
and line defects (sourced by E^jSLJ 



B. The bF + bb description of the 3D semion model 

Having seen that (3+l)-D bF theory gives an effec- 
tive description of topological order in the 3D toric code, 
we now seek to motivate a field theory which has been 
claimecPto give a similar effective description of 3DSem. 
Concentrating on 3DSem on the solid donut (see Sec. IV I, 
we expect the field theory to capture two important fea- 
tures of topological order: the ground state degeneracy 
(which was 2), and the associated deconfined semionic 
excitations on the surface of the manifold (which is 
T 2 ). These are precisely the topological properties of 
a v = 1/2 bosonic Laughlin state on T 2 . Moreover, it 
is well known that the topological properties of v = 1/2 
Laughlin state are described by a £7(1)2 Chern-Simons 
theory. Therefore, in the absence of sources in the bulk, 
we expect our effective (3+l)D field theory on the solid 
donut to reduce to a (2+l)D ?7(1)2 Chern-Simons theory 
on the boundary T 2 . 

In fact, in the lattice model the surface string opera- 
tors (W(C Z ) and W(C$)) used to toggle and differenti- 
ate between the different ground state sectors have exact 
analogues in the field theory. Because the semion graph- 
ical rules Fig. [2 are based on a U(l) 2 abelian Chern- 
Simons theory, the non- contractible string-operators on 
the surface obey a similar algebra to the non-contractible 
Wilson-line operators that determine the dimension of 
the ground state Hilbert space in an abelian Chern- 
Simons theory living on T 2 . 

This suggests a (3+l)D field theory with action 



& v (x) =J d 2 a5^(x-X{<j)) 



dX p dX v dX p dX v 



da 2 



da 2 da 1 

where {X p (a 1 , a 2 )} is an embedding of a 2D world-sheet 
into the 4D space, parameterized by (<r , a 2 ). The point 
particles and line-like objects of the field theory corre- 
spond respectively to the vertex defects and lines of pla- 
quette defects in the lattice model. 

We expect the field theory to give us the ground state 
degeneracy of the 3D toric code, as ground state degen- 
eracy is an important topological feature of the lattice 
model. More precisely, the dimension of the Hilbert space 
of the field theory in the absence of sources should equal 
the ground state degeneracy of the lattice model. Indeed, 
it can be shown that the Hilbert space corresponding to 



the bF action in (21 1 has dimension 2 3 on the 3-toru92B, 
in accordance with the lattice calculations in Sec. lIII ATI 



Sff [A] 



d x f c p ^ F n V F pa 

\ l07T 



, (24) 



with k = 2 and Fp V = dpA v — d v Ap, because it reduces 
to a surface Chern-Simons term using e pL /p ° 'F \ IV F y pa — 
4 dp (e' lvpo A v d p A a ) to give action 



Sff [A] = — I ,11 



dM 



dSnp^ vpa A v d p A a + \ d 4 xApJ p , 



ScslA] 

(25) 

where dS is the area element on dAA , and Up is the outer 
normal. The first part of this expression is just a Chern- 
Simons action of the gauge field on the surface of the 
manifold. Heuristically, the total partition function takes 
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the form 
Zff [J] 

= J DA s e^ Scs ^ +i f Js - As J DA B ^ Jb ' Ab 

= J DA s e^ Scs[As]+ ^ ^-^S^ = in bulk] , (26) 

where we have smoothly split A^ and J M into bulk (B) 
parts and surface (S) parts. J M can serve as a source term 
for semionic excitations on the surface. However, when 
one tries source a particle world-line in the bulk using 
a J M term the field integral disappears upon integrating 
out the bulk degrees of freedom of A^ , as represented by 
the delta function constraint. 

This suggests that the Sff action gives an incomplete 
description of 3DSem, because we saw in Sec. HI B 2| that 
point defects do appear in the bulk and moreover they 
are attached to the end-points of line defects. It would 
be reassuring to see these compound 'point + line' ob- 
jects emerge naturally from a field theory. To this end 
we introduce a source E^ and a corresponding rank-2 
antisymmetric field &„„ to give a 'bF + bV action 



S bF+bb [A,b] = J fad""" (A.b^dpA, 



I67r b ^ bpa 



(27) 



where J p sources point particle world- lines (Eq. (22 1), 



and E^ sources the world-sheets of line defects 
(Eq. p3j)). Upon integrating out b^ we are left with 
an effective action 



SbF+bb [A] - ^Scs [A] 



d 4 x 



e^E^E^ + A^^ + d^) , (28) 



where Scs is again the Chern-Simons action on the sur- 
face of the manifold and we have assumed for now that 
E vanishes on this surface. In the absence of E sources 
this reduces to the Sff action in Eq. (25), therefore the 
Hilbert space of the field theory will be isomorphic to 
the ground state subspace of the 3D semion model. Fur- 
thermore, J can again be used to source semions on the 
surface of the manifold. Leaving the £ sources in, we can 
further deduce (heuristically) that 



Z [J, S] 



e 4k 



us; 



DA s e* 



:Scs[As]+if Js-As 



DA B e JV+*M* 



/ DA £Scs[As]+ifJs-As S rj = 



-if in bulk] 
(29) 



where K» = d u E» v , and I {E) = / d^e^E^E^. In 
the equation above we have smoothly split A into a bulk cal rules, taking the general form shown in Fig. 21 we 



(B) part and a surface (S) part, the only restriction on 
this splitting is Ab — on the boundary. The delta 
function constraint in Eq. (29) ties bulk vertex defects 
(J^) to . It is a general feature (see example below) 
of line sources that — everywhere but at the ends of 
the line source. With this observation, the delta function 
constraint forces bulk vertex defects (J M ) to lie at the 
end-points of line sources (where ^ 0). 

As an example, let us look at a line defect in the z- 
direction which extends from (0, 0, —a) to (0, 0, a). Using 
the prescription Eq. ( 23 1 , it can be shown that the E^ v 
source takes the form 

E"" = I(z G [-a, a]) S (2) {x, y) (5^ - S^S*) (30) 

K" = <5 (2) Or, y) (5(z + a) - S(z - a))^ t (31) 

where I[z G [—a, a]) = 1 when z G [—a, a] and I(z G 
[—a, a]) = otherwise. Note that is zero everywhere 
but at the endpoints of the line source. Therefore, from 
the delta function constraint in Eq. ( [29] ), point charges 
are found precisely at the end-points of the line source 
i.e. = ( x , y ) (5{z -a)- 5{z + a))8^ t . 

Therefore bulk point particles are found precisely at 
the endpoints of line-like defects where ^ 0; the 
only other place we can have non-vanishing J M is on the 
boundary of the manifold. This corresponds to what we 
found on the lattice in Sec. |iVB| point defects in the bulk 
always lie at the ends of line defects of plaquette viola- 
tions because point defects are confined, but point defects 
on the surface of the manifold need not lie at the ends 
of line defects. We have shown that the bF + bb and FF 
field theories capture important aspects of the topological 
order in the lattice model, namely the ground state de- 
generacy and the surface anyonic excitations. However, 
we also showed that, unlike the FF theory, the bF + bb 
field theory can be used to source the bulk 'line+point' 
defects characteristic of the 3DSem lattice model. 

While the FF and bF + bb field theories can describe 
the topological features of the lattice model, they also 
contain extra gapless degrees of freedom not present in 
the lattice model. These extra degrees of freedom are 
not apparent in either field theory as both have vanish- 
ing Hamiltonians; however, introducing a Maxwell (or 
Yang-Mills) term in the bulk will render these degrees of 
freedom dynamical, gapping many of these formerly gap- 
less degrees of freedom but leaving a gapless photon in 
the bulk. In order to obtain a field theoretic description 
that faithfully captures the low-energy (gapped in the 
bulk physics of the lattice model, these gapless degrees 
of freedom must be eliminated from the field theory22. 



VI. GENERAL WALKER- WANG MODELS 

So far we have studied two of the simplest Walker- 
Wang models, built from the graphical rules in Fig. [2] 
There are many more examples of Walker- Wang models, 
each based on a different set of self-consistent graphi- 
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refer to such a collection of rules as a unitary braided 
fusion category (or just 'category' for short). In this sec- 
tion we introduce the general Walker- Wang models in 
Sec. |VI A| and find that, of the many possible models, 
there is a family which behaves in a similar manner to 
3DSem. We call these 'confined Walker- Wang models' 
because, like 3DSem, all of their bulk excitations are con- 
fined. In Sec. |VIB| we show that a Walker- Wang model 
is 'confined' if and only if it is based on a modular tensor 



(a) Vertex 
constraints 



category (MTC) (which we define in Sec. VI B I. 



We then go on to demonstrate that confined Walker- 
Wang models share many other properties with 3DSem. 
First we show that the confined models have a unique 
ground state on the 3-torus in Sec. |VIC"T) We then show 
in Sec. |VI C 2| that confined models can have multiple 
ground states on manifolds with boundary like the solid 
donut, and these degenerate ground states are associ- 
ated with deconfined anyonic vertex excitations which 
are pinned to the boundary. These results support the 
suggestion of Walker and Wang^ that the lattice models 
based on certain MTC's (specifically those MTC's com- 
ing from a Chern-Simons theory) have a description in 
terms of bF + bb theory, analogous to the bF + bb descrip- 
tion of 3DSem in Sec. IVBl 



The general Walker- Wang Hamiltonian, and 
excitations 



The Hilbert spaces of the 3D toric code and 3D semion 
model consisted of a 2 state system <j z = ±1 on each 
edge of the lattice shown in Fig. [lOj and we represented 
spin configurations graphically by coloring in only the 
<j z = — 1 edges. The ground states were determined 
by using the graphical rules: Fig. |2ja) told us that the 
ground state involved configurations with an even num- 
ber of colored edges coming into a vertex (so the ground 
state was a 'loop gas'), while the graphical rules Fig.pFb)- 
(e) provided relationships between the amplitudes of dif- 
ferent loop gas configurations. 

More general Walker-Wang models are based, in an 
analogous way, on more complicated sets of graphical 
rules. For example, probably the next simplest Walker- 
Wang model is based on the Fibonacci category. The 
graphical rules for this model are given in Fig. [20](a)-(e). 
The Walker- Wang model based on this category(3DFib) 
still consist of a 2 state system <j z — ±1 on every edge, 
and again we graphically represent different spin config- 
urations by only coloring in edges with <j z = —1. The 
Hamiltonian (which we define below) is set up so as to 
ensure the following: The ground state is a superposition 
of configurations obeying the modified vertex constraints 
a), and the graphical rules Fig. 



Fig. |20[a), and the graphical rules Fig. |20[b)-(e) give 
new linear relations between the amplitudes of different 
kets in the ground state superposition. As for the semion 
rules, the Fibonacci rules Fig. [20^a) permit vertices with 
none or two of their incoming edges colored, and for- 
bid vertices with a single colored edge. However, unlike 




(f) u>q - strin 



FIG. 20. The figure above shows the graphical rules for the 
Fibonacci category, where 7 = (1+ v / 5)/2 and £ = e l7T ^ 5 . (a) 
represents the fact that the ground state of a Walker-Wang 
model based on this category is a superposition of configura- 
tions involving vertices with either zero, two or three edges 
colored in. The rules (b)-(e) give graphical relations between 
the amplitudes of different spin configurations on the ground 
state, as well as provide a graphical calculus for defining string 
operators. Note that there is a rule conjugate to (e) obtained 
by turning the over-crossing into an under-crossing on the left 
hand side, and sending £ — > f* = £~ on the right hand side. 
The string type ujq in (f) is used in the definitions of plaquette 
projection operators. 



the semion rules, the Fibonacci rules also permit vertices 
with all three incoming edges colored. This means that 
the ground state does not quite look like a loop gas - 
we say that the ground state is a 'string-net^ obeying 
Fig.^a). 

Further generalizations of Walker Wang models can in- 
clude Hilbert spaces where each edge is in one of k pos- 
sible states, as compared to the above cases where each 
edge is in only one of two possible states. More generally, 
a Walker Wang model can be defined by any 'unitary 
braided fusion categories^. A 'unitary braided fusion 
category' is a set of possible edge labels along with a set 
of graphical rules (sh own in Fig. 21 1 that satisfy certain 
consistency condition d 25 * 26 ! which we will not enumerate 
here. For a category with k possible labels, the Hilbert 
space of the corresponding WW model consists of giving 
a label i G {0, . . . , k — 1} to each edge in the system^. 

The Hamiltonian, which we introduce in the next sec- 
tion, is constructed to ensure two things. Firstly it en- 
sures that the ground state is a string-net with k — 1 dif- 
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ferent 'colored' strings with labels {1, . . . , k — 1} obeying 
a more general set of vertex constraints represented by 
Fig- 
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'&). Secondly, it ensures that the rules Fig. 21 'b)- 
(e) give linear relations between the amplitudes of differ- 
ent string-net configurations in the ground state super- 
position. 



1. Hamiltonians 

Having explained that all Walker- Wang Hamiltonians 
are contrived so as to ensure the ground state is described 
by a set of graphical rules, we now define the Hamiltoni- 
ans. The details in this section are not essential to un- 
derstanding the qualitative behavior of confined Walker- 
Wang models, which we resume studying in Sec. |VI~C| 
The casual reader may therefore wish to skip directly to 
Sec. IVICl 

The general Walker- Wang Hamiltonian takes the form: 



H 



WW 



(32) 



where V V ,V P are a commuting set of projectors. The 
ground state(s) can be solved for exactly, and has the 
defining property that V v — V p — 1 for all vertices and 
plaquettes. The operator V v projects onto configurations 
where only certain combinations of edge colors come into 
a vertex (see Fig. [21~|a)), and so the condition V v = 1 
for all v ensures that the ground state is a string-net 
involving only these types of configurations. The plaque- 
tte operators will be defined below, but we reiterate that 
they are constructed to ensure that Fig. [2l|b)-(e) relate 
the amplitudes of different string-net configurations. 

All that remains is to define the plaquette projectors 
generally. For explicit expressions for the V p , we re- 
fer readers to [8]; in what follows, we will define the 
V p using a graphical calculus. We saw in the 3D toric 
code and 3DSem models that there exists a graphi- 
cal mnemonic for acting with string-operators. Ex- 
actly the same mnemonic works here. W V (C) denotes 
an operator formed by laying down a string of type 
i G {0, 1, . . . , k — 1} along a path C and fusing it into 
the edges using Fig.[2l|b)-(e). While there are only fc — 1 
colored strings in the category, it is notationally conve- 
nient to introduce an 'empty' string labelled i = with 
the convention that it braids and fuses trivially with ev- 
erything, and has A^ = o = 1 (see Fig.[2l|c)). This implies 
that empty string operators always acts trivially on the 
Hilbert space i.e. W v =0 (dp) = I. 

We now apply this formalism to define the plaquette 
operators. First we define plaquette operators W v (dp) 
for each string type i: Using the fixed projection of the 
lattice onto 2D, thread a string carrying label i around 
dp using the under and over-crossing prescription shown 
in Fig. [l4|(d)-(f) and fuse it into the edges of p using the 
graphical rules (Fig. 21 'b)-(e)); this prescription ensures 



that adjacent plaquette operators commute, as shown in 
Fig. [22j Now, the plaquette projector is defined by 



1 



fc-i 



j=0 



(33) 



where the A 4 is known as the quantum dimension of i 
(see Fig. [21 c)), and I? 2 = X^~o A f ■ This operator will 
sometimes be referred to as Wy°(dp), where ujq is the 
superposition of string types 



w 



fc-l 

E 

i=0 



Xi/v 2 



(34) 



where i indicates a string of type i. While the form of 
this ujq operator may seem complicated, it is constructed 
to have one crucial property namely the 'handle-slide' 
property 2 ^ shown and explained in Fig. 23 It is precisely 
this property which ensures that Fig. 21 b)-(e) relate the 
amplitudes of different configurations in the ground state. 



2. DSem and SDSem in the more general language 

Let us now try to understand 3DSem in this new lan- 
guage. The 3DSem category has two string types. There 
is the 'empty' string carrying label i = 0, which for sim- 
plicity we did not mention in Fig. [2] because it braids 
and fuses trivially with everything. Then there is the 
single colored string type which has label i = 1 and has 
semionic self-braiding Fig. |2je) . 

The allowed vertices in 3DSem had only an even num- 
ber of down spins attached to each vertex, which corre- 
sponded with B v = +1, where we defined B v in Eq. (14 1. 
Therefore we should define V v = (1 + B v )/2 so that the 
condition V v = 1 is equivalent to B v = 1, a condition 



which ensured that the ground state was a loop gas (see 
Sec. |HTB ). 

According to Eq. ( [33] ) the definition of V p in DSem is 
V p = (1 - W v {dp))/2, because A = 1, Ai = -1. But 



we noted in Sec. Ill B that B p has a graphical definition 



which is precisely the definition of the string operator 



W v (dp). Hence V p = (1 



B p )/2. Plugging these rela- 
f B v )/2, into Eq. ([32]) gives 



tions, as well as V v = (1 
back Eq. ( 14 1 up to an unimportant re-scaling of the 



Hamiltonian (and an overall shift of the Hamiltonian by 
a constant). 



B. Modular tensor categories and excitations 



In Sec. lIIIB~2l we found that the bulk vertex excitations 
of the 3DSem model were confined, and this was related 
to the statement that string operators have an energy 
cost for every plaquette they thread. For these reasons, 
we called 3DSem a 'confined Walker- Wang model'. We 
now turn our attention to Walker- Wang models based on 
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General unitary braided fusion category 


(a) Vertex 
constraints 




C 1, if (a, fc, cj allowed 
c = \ 0, otherwise. 


(b) Deformation 




a ^J^a 


(c) Loop 
collapsing 


0° 


(d) Fusion 


W - iz>th "Y 

a' \d 


(e) Braiding 
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r . , Y 


(f) too - string 




(g) S- matrix 





FIG. 21. (a) Represents the vertices allowed by the cate- 
gory; the ground state of a Walker- Wang model will involve 
only these types of vertices. The diagrams in (b)-(d) serve 
two purposes. Firstly, they tell us the relative amplitudes of 
'string-net' configurations in the ground state e.g. row (c) 
tells us that configurations related by removing a closed loop 
carrying label a occur with a relative factor of A a in the 
ground state. Second, these diagrams provide a neat graph- 
ical mnemonic for the definitions of string operators. Note 
that there is a rule conjugate to (e) obtained by turning the 
over-crossing into an under-crossing on the left hand side, and 
sending R — > R* on the right hand side. 



modular tensor categories (MTC's). We will show that 
these are precisely the Walker- Wang models which are 
confined, in the sense that their string operators have an 
energy cost for every plaquette they thread. 

MTC's are categories with a unitary 'S- matrix^. 
They have the special property that a loop of the 
string type used to define V p (i.e. ujq = ^ Ji ^ii/'D 2 in 
Fig. [2r}f)) kills all flux going through it i.e using the fu- 
sion and braiding rules in Fig. 21 to evaluate Fig. [23^b) 
yields an answer of zero if the flux i ^ 0. Examples 
of MTC's include the semion, Fibonacci, and Ising cat- 
egories as well as SU(2) for all rp3. An important 
non-example is the toric code; we see in Fig. 23 [c) that 
a string operator (1 + W^)/2 does not kill flux going 
through it. 



an i string along C, and fusing it into the edges along the 
path. This produces vertex defects at the end-points of 
C. For the 3D toric code, such string operators only pro- 
duced vertex defects at their endpoints. But for 3DSem, 
the operator also produced plaquette violations along 
C, which represented the fact that bulk vertex excita- 
tions are confined in 3DSem. We will now show that 
MTC based models behave like 3DSem: Acting on the 
ground state of any MTC based model with string op- 
erator Wy(C) will excite plaquettes p threaded by the 
string in addition to creating vertex defects at the ends 
of C. Put another way, point defects are confined because 
they are tied to a line of plaquette defects. 

The previous paragraph establishes that vertex string 
operators that act in accordance with the graphical rules 
produce confined excitations. In Appendix [B] we show 
that this holds for any string operator in a confined 
Walker- Wang model, completing the proof that all ex- 
citations are confined in the bulk. 

Algebraically, confinement means that V p Wy(C) 
GS) = for i =/= 0, where p is any plaquette threaded by 
C. We prove this relation graphically in Fig. [24^ a)- (d). 
The crucial feature of MTC's is that loops of too string kill 
the flux going through them (see Fig. 23 1 , and this key 
property is used to evaluate Fig. 



24|d)7^ This behavior 



should be contrasted with that seen for the toric code, 
which is not a modular tensor category. For the toric 
code we can follow the same reasoning as for a MTC in 
Fig. 24'a)-(d). The difference comes in evaluating dia- 



gram"^) ; in the toric code the ojq string does not kill the 
i flux, it gives back the same answer Fig. 24 'a) because of 
the identity shown in Fig. [23|(c) . Algebraically this can 
be expressed as V p W l | GS) = W l | GS), and this im- 
plies that the toric code string operator does not produce 
any plaquette violations along its length, a point we have 



We can define a string operator W V (C) by threading 



already discussed in a simpler setting in Sec. |IIIB 2 

Having seen that all MTC based models are confined, 
we ask the reverse question: If a category is bulk con- 
fined, is it necessarily modular? The answer is yes. If 
a category is not modular, then it can be showrP3 that 
there must exist a non- vacuum particle which braids triv- 
ially with all other particles. If we write down a string 
operator carrying the label of such a particle, we will 
find that it commutes with every plaquette projector, 
and therefore gives rise to deconfincd bulk point excita- 
tions. Therefore 'not-MTC implies 'not bulk confined', 
or equivalently all bulk confined models must be based on 
MTC's. With the reasoning in the previous paragraph, 
this implies that a model is bulk confined if and only if 
it is based on a MTC. 



Confined Walker- Wang models on various 
manifolds 



We have shown that string operators in confined 
Walker- Wang models violate every plaquette they thread. 
This suggests that confined Walker-Wang models on 
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FIG. 22. This figure shows the graphical definition of plaque- 
tte operators in general Walker Wang models. Having fixed 
a projection we define a plaquette projector by threading a 
ujq string just inside the perimeter of the plaquette. Two ex- 
amples are shown in (a) and (b). Then in (c) we show that 
these adjacent operators commute with each other because 
the diagrams representing V pl Vp 2 and V P2 'Pp 1 are identical. 



manifolds without boundary have a unique ground state, 
an assertion we verify in Sec. |VI C 1| (this result is known 
to the mathematical community 421 ^ . We then go on to 
show in Sec. 



VIC 2 that, like 3DScm, confined Walker- 



Wang models on manifolds with boundary have decon- 
fined anyonic excitations pinned to the boundaries, and 
that these deconfined excitations are associated with a 
ground state degeneracy. 



The 3-torus 



One consequence of bulk confinement is that there 
are no non-contractible string operators in the bulk 
that commute with the Hamiltonian of an MTC based 
model. This naively suggests that the ground state is 
non-degenerate on manifolds without boundary because 
there are no operators which commute with the Hamil- 
tonian and toggle between different ground states. We 
have already written out an argument specific to the 3D 
semion model, proving that it has non-degenerate ground 
states on the 3-torus. In Appendix |A"| we present a gen- 
eral proof that the ground state is non-degenerate on T 3 
for any modular tensor category. For MTC Walker- Wang 
models it appears as though the interesting physics oc- 
curs on the surface, and not in the bulk. This motivates 
the study of a simple manifold with boundary. 




r- 




Modular tensor category 



The Toric Code 



(b) 



(c) 
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1 
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FIG. 23. (Color online) (a) demonstrates the handle-slide 
property of the ujq loop which appears in the definition of 
plaquette operator. A special property of cjo allows us to 
deform other string operators (or indeed edge strings) i over 
ljo regardless of whether other strings thread luq. The X in 
the middle of the loo loop indicates that other strings may 
be passing through the middle of the loop, (b) In a modular 
category no string (except the empty string) can pass through 
an luo loop. However, in a non-modular category, such as in 
the toric code, this is not the case. 



2. The solid donut D 2 x S 1 

We now put a MTC based model on a solid donut 
choosing the 'smooth' boundary conditions shown in 
Fig. [25j These are completely analogous to the smooth 
boundary conditions described for the 3D toric code and 



3DSem model in Sec. IV and are similarly constructed 
to gap the surface to vertex and plaquette excitations. 

As in the case of 3DSem on the donut, we can create 
pairs of deconfined anyons on the surface of the mani- 
fold by using graphically defined operators of the form 
Wy (Cab), where Cab is a path lying just above the 
surface of the manifold connecting points A to B and 
i € {1, . . . , k — 1} is any label in the category. The 
anyons are deconfined because the path of the string does 
not thread any plaquettes (see Fig. 25 'd)). 

For such an MTC with k — 1 non-trivial string 
types, the ground state degeneracy for the correspond- 
ing Walker- Wang model on the solid donut is k - exactly 
that of the chiral anyon theory living on its surface. To 
see this, note that any ket without vertex violations can 
be related, using graphical rules Fig. [21] to a combina- 
tion of the k canonical kets shown in Fig. 25 ^b). The 



ground states labelled by these canonical kets are mu- 
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Toric code: V p W P p | GS) =W P p | GS) 
MTC: 

PpWPp \ GS) = 8 iWPp |GS) 

FIG. 24. (Color online) This diagram shows that acting on 
the ground state of a lattice model with a string operator 
carrying label i 7^ and threading plaquette p will excite the 
plaquette to V p = if the underlying category is modular, 
but will not excite the plaquette if the underlying category 
is Z2 (the toric code case), (a) Represents the ground state, 
which we have acted on with V v which equals 1 on the ground 
state. In (b) we act with a string operator. In (b) is we act 
with plaquette projector V v to ascertain the new eigenvalue 
of V p and this forms (c). Using the handle-slide property 
of the plaquette projectors strings, we handle-slide the inner 
projector over the outer projector. The result is (d) in which 
a c^o-string is left encircling the i-string giving Sot for a MTC, 
which is zero because of our choice of i. 



tually orthogonal; we can verify this by showing that 
they have distinct eigenvalues under flux measurement, 
which involved acting with a close d s tring operator of 
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d). Hence there 



the form Wy l (C^) defined in Fig. 
are precisely k distinct ground states, and we can toggle 
between them usi ng closed string operators of the form 
Wy (C z ) (see Fig. 25 'c)). This is consistent with the de- 
generacy of 2 founa for 3DSem in Sec. |IV| As a further 
example, an MTC based on SU(2) n has a degeneracy of 
n + 1. We note again that the ground state degeneracy 
is tied to the existence of surface anyon excitations, a 
point which motivates the field theoretic description in 
the next section. 



D. The bF + bb description of MTC's 



We saw in Sec. VB that a level k = 2 U(l) bF + bb 





theory captured the ground state degeneracy and types 
of defects of the 3DSem model. We now turn our atten- 



FIG. 25. (a) Shows the k independent ground states | GS%) 
where i £ {0, 1, 2, . . . , k — 1} and V is the ground state projec- 
tor. The dotted lines are edges in state while the bold black 
edges are in states corresponding to the i-label. (b) Defines a 
string operator W v (C z ) winding around the non-contractible 
2-cycle of the solid donut. This operator toggles can be used 
to toggle between ground states e.g. W v (C z ) | GSo) =\ GSi). 
(c) Defines yet another operator Wy z (C^) which can be used 
to distinguish between the different ground states. This op- 
erator involves a string type Wj with a special superposition 
of labels Ui = ~^2 h S^h where S is the S-matrix of the un- 
derlying category. In a MTC, this operator has the property 
Wy*- (Cj,) I GSj) — Sij I GSj). (d) Defines an open string op- 
erator, creating defects at A and B but otherwise commuting 
with the Hamiltonian along its length. 



tion to the more general statement If a Walker- Wang 
model is based on a MTC which looks like a Chern- 
Simons anyon theory, for example SU(2)„™ then the 
resulting lattice model has an effective description as a 
bF + bb field theory with gauge group SU(2) and level n. 

For concreteness, let us fix our attention on the mod- 
ular Walker- Wang models based on SU(2) n . From our 
work in Sec. |VI C\ we see that these models have a ground 
state degeneracy of n 4-1 on the solid donut, which is pre- 
cisely the dimension of the Hilbert space of a non-abelian 
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SU(2) n Chern-Simons theory. Furthermore, these mod- 
els have deconfined anyons on their surface, which obey 
the statistics found in the SU(2) n Chern-Simons theory. 
This information so far suggests the non-abelian effective 
field theory 



Sff [A] = J d*x Tr (JLe^F^Fpo') , 



(35) 



which reduces to the correct surface Chern-Simons the- 
ory. Above we have assumed that A lives in the fun- 
damental representation of the gauge group, and F is 
the corresponding non-abelian field strength. Creat- 
ing surface nonabelions corresponds to introducing path- 
orderecpSl Wilson lines of the form 



(36) 



where 1Z is a representation of the gauge group, which 
labels the particle type. However, using the same heuris- 
tic reasoning that led to the <5-function constraint in 
Eq. (261, this field theory seems to forbid the presence of 
bulk point particles. As we saw in Sec. |VI B[ bulk point 
particles do exist in the lattice model, but are tied to the 
ends of lines of plaquette defects. This again motivates 
the introduction of a two form b^ transforming in the 
adjoint representation to give 



S bF+bb [A, b]=J d 4 xne> iU <" T Tr (J^b^F^ - -j^M 



pa 



(37) 



which reduces to Eq. ( 35 ) upon integration over b, so re- 



produces the correct ground state degeneracy and surface 
anyonic statistics. To source a compound 'line+point' de- 
fect, one introduces Wilson line and surface operators of 
the form 

Tr K , Lif^^ Tr K ( e * f (38) 

There are many subtleties in defining surface 
operators®! which we do not discuss here. That aside, 
we expect that for such a 'line+point' defect to exist in 
the bulk, the following analogue of the delta function in 
Eq. p9h arises: 8 n >n5 [J = -K = -d v W v in bulk]. In 



other words, the point defects are bound to the ends of 
line defects, and both carry the same quantum label. 



VII. CONNECTION TO MODELS WITH 
SURFACE FRACTIONAL QUANTUM HALL 
STATES 



Laughlin states. In particular we concentrate on a model 
which, on its surface, has topological order resembling a 
v = 1/3 fermionic fractional quantum Hall effect. 

Though the microscopic degrees of freedom of Walker- 
Wang models are fc-state spins (i.e., k possible quantum 
numbers on each edge), their low-energy properties can 
be made to resemble those of fermionic systems. The 
resemblance that we will describe is reminiscent of that 
between the 2D toric code and a superconductor in two 
dimensions^. Though the Toric code does not describe 
superconductivity per se, both systems have a ground 
state degeneracy of 4 on the torus, and two types of low- 
lying excitations (vertex and plaquette for the Toric code, 
or BdG quasi-particles and ir vortices for the supercon- 
ductor) which acquire a Berry phase of — 1 when braided 
around one another. In other words, the 2D thin-film 
superconductor and the 2D toric code exhibit the same 
kind of topological order^U 

In this section we give examples of lattice models whose 
surfaces exhibit the same kind of T-symmetry breaking 
topological order as a fermionic fractional quantum Hall 
state. The existence of these phases is suggestive: In 
the same way that the toric code has the topological or- 
der of a superconductor, the lattice model in Sec. |VIIB| 
may well have the same topological order as some sys- 
tem of electrons. We will not discuss such electron sys- 
tems in detail, but we note that the lattice model in 
Sec. |VII B| shares some of the topological properties one 
might expect of a fractional topological insulator under- 
going confinement^, namely, time-reversal symmetry is 
explicitly broken, leading to surface states with the topo- 
logical order of a fermionic Laughlin state, and the only 
bulk deconfined excitation is the fermion. 

This section is structured as follows. In Sec. I VII Al we 
revisit the 3DSem model, which has a confined bulk and 
the topological order of a v = 1/2 bosonic fractional Hall 
effect on its surface. Then in Sec. |VII B| we attempt to 
find a lattice model which has surface topological order 
resembling a fermionic v = 1/3 fractional Hall effect. 
We succeed in part: we find a lattice model with the 
expected deconfined anyons pinned to its surface, but we 
also find that there exists a deconfined vertex defect in 
the bulk with fermionic statistics. In Sec. |VIIB 5| we give 
a dictionary for going between more general fermionic 
fractional Hall effect, and a Walker-Wang model with 
the corresponding topological order on its surface. Once 
again we find that lattice models with a fermionic surface 
fractional Hall effect have a deconfined fermionic vertex 
excitation in the bulk. 



A. Bosonic v = 1/2 surface Hall state 



We have seen that the 3DSem model has a surface that 
is topologically a bosonic Laughlin state at 1/2-filling (i.e. 
described by a k — 2 abelian Chern-Simons theory^, as 



discussed in Sec. VB|. In this section we present lattice 



models with the surface topological order of fermionic 



We saw in Sec. IIVB II that the 3D semion model has 
surface excitations with semionic statistics - in other 
words, with exactly the statistics of the charge-1/2 ex- 
citations of the bosonic v = 1/2 Laughlin state (which 
are described by a U(l) Chern-Simons theory at level 



2G 



k = 2jPH. Because excitations in the bulk are confined 
and the ground-state degeneracy is determined only by 
the surface, the topological order of 3DSem on the (3D) 
solid donut is exactly that of the bosonic v — 1/2 Laugh- 
lin state on the torus. 

We begin by making more precise the analogy between 
the surface of the lattice model and a system of bosons 
of charge qs in a magnetic field. The charge gs/2 excita- 
tions have the same semionic statistics as surface vertex 
defects (corresponding to vertices at the end s of the string 
operators W(Cab) discussed in Sec. IV B 1 ). The bosonic 



charge qs excitations, however, do not exist in the lattice 
model: excitations that carry charge but no interesting 
statistics are, in the string-net framework, indistinguish- 
able from the identity. This represents an important dif- 
ference between the Walker- Wang lattice models of this 
section and actual quantum Hall fluids: The lattice mod- 
els have Z n rather than U(l) charge conservation. In the 
present example fusing together two of the <z_b/2 exci- 
tations in the bosonic v = 1/2 state gives a charge qs 
excitation, which is physically measurable but topolog- 
ically trivial, while in 3DScm fusing the corresponding 
vertex defects gives the identity - so the lattice model 
has only Z2 charge conservation. 



B. Fermionic v = 1/3 surface Hall state 

We have seen that the 3D semion model has the same 
topological order as a v = 1/2 bosonic fractional Hall 
effect on the boundary of a confining bulk. Can we con- 
struct a Walker- Wang model with surface topological or- 
der resembling a v = 1/3 fermionic Hall effect? We will 
find that there does exist such model but, unlike 3DScm, 
it has a deconfined vertex excitation (with fermionic self- 
braiding) in the bulk which is associated with an addi- 
tional ground state degeneracy. Therefore, the lattice 
model does not quite have the topological order of a sur- 
face v = 1/3 Laughlin state on a confining bulk. This 
difference can, however, be eliminated by adding a term 
to the bulk Hamiltonian. 

We choose a category with labels which correspond 
with the quasi-particle content of the v — 1/3 quan- 
tum Hall effect (summarized in Table |l| . The category 



is known as Zg 2 ' 1 x Z^, and it can be thought of as aris- 
ing from a compact k = 3 U(l) Chern-Simons theory^. 

The category in question has labels represented by 
a doublet i — (11,12) where i\ is an integer modulo 3 
and 12 is an integer modulo 2. The Hilbert space thus 
consists of a 6-state system on each edge of a lattice^. 
The label (ii = 1,22 = 1) represents a Laughlin quasi- 
particle of charge e/3, while (ii = 0,i2 = 1) represents 
a charge e quasi-particle (which, unlike the boson of the 
v = 1/2 Laughlin state, we must keep track of explic- 
itly because of its fermionic statistics). Excitations such 
as (ii = l,i2 = 0) can be viewed as the combination of 
Laughlin charge e/3 and charge e quasi-particles. The 
correspondence between the doublets and the various 
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FIG. 26. (a) Shows the six possible edge labels of the cat- 
egory. All of the string labels should be viewed as coming 
from x Zj 1 '. (b) to (f) show the category rules, which 
we have explained in Sec. |VI| Note that labels can be added 
component by component e.g. a + b — (ai +61 mod 3, 02 + ^2 
mod 2). (g) Shows the Wo-string in the category (which, as 
always, has the handle-slide property), while (h) shows the 
S-matrix of the theory (where r\ = e* 2 "^ 3 ). The columns and 
rows of the S-matrix are ordered according to the list in (a). 



Laughlin quasi-particles is summarized in Table [T] 



Excitations of the Zg 2 ' x Zj 1 ' model 



In this section we summarize the forms of excitations 
present in the Walker- Wang model based on Z3 x Z 2 ^. 
First we discuss the vertex defects , which carry a Z 6 = 



Z^ 2) x Z 2 1; charge (see Sec. IVIIAJ). Unlike the confined 
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Walker- Wang models such as 3DSem, the Walker- Wang 

has deconfined point-like ex- 



model based on Z 



(2) 



7(D 

J 2 



citations in the bulk and on its surface. The deconfined 
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Charge 


Exchange Phase 


Category Label 





1 


(0,0) 


e/3 




(1,1) 


2e/3 


e i4.w/3 


(2,0) 


e 


e 17 * = -1 


(0,1) 


4e/3 


e i4.w/3 


(1,0) 


5e/3 


e in/3 


(2,1) 



TABLE I. This table shows the edge labels in the category 
x Zj 1 ' used to construct the state with surface resembling 
a v — 1/3 fermionic Hall effect. Each of the labels in the 
category correspond with a quasi-particle in the Laughlin v = 
1/3 state. The exchange phase is the phase accumulated by 
exchanging two particles of the same type. 



bulk excitations are fermionic vertex defects that occur 
at the end-points of strings labelled (0, 1). The other four 
possible vertex defects, which on the surface correspond 
to Laughlin quasi-particles, or Laughlin quasi-particles 
bound to electrons, are deconfined only on the surface. 
This is consistent with the fact that point particles in 3D 
must be bosons or fermions. 

It is enlightening to see how this arises in practice. 
From the S'-matrix (Fig. [26^h)) we see that strings car- 
rying label / = (0, 1) (/ for fermion) braid trivially with 
all other particles (represented by the fact that the col- 
umn Sif is identical to Sio). This in turn implies that 
string operators of the form W V (C) commute with all 
the plaquettes along their length because they commute 
with all the string operators used to define the plaque- 
tte operator. On the other hand, all other string types 
braid non-trivially, so there is only one deconfined vertex 
defect. Thus the other four types of excitation listed in 
Table [I] are confined in the bulk, but as in the case of 
MTC based models, they are deconfined on the surface 
of the manifold. 

In the bulk, the four kinds of confined vertex defect 
are ass ociated with lines of plaquette defects (like those 
in Sec. 



IIIB2) which we can think of as 'chiral vortex' 



lines. However there is another type of confined vortex 
line which has the same relative statistic with the de- 
confined (0, 1) vertex defects as we found in Sec. Ill A 
between vortex lines and vertex defects in the 3D toric 
code. We can define this vortex line using a surface op- 
erator (see Eq. (13)) of the form 



W P {S) 



(39) 



where <rj acts on edge j in the following way: if edge 
j is in state i = (11,12) then aj = ( — l)' 2 . Here S is a 
surface on the dual lattice, and the product is over edges 
cutting the surface transversally. The operator Wp(S) 
creates a line of plaquette violations on its boundary in 
the manner shown for the 3D toric code in Fig. |13| In 



3D toric code: the deconfined (0, 1) vertex defects behave 
like the vertex defects in the toric code (except (0, 1) have 
fermionic self-braiding), and the vortex lines created by 
Eq. ( 39 1 are like the vortex lines in the toric code. Having 



understood the spectrum of this Walker- Wang model, we 
now investigate its properties on some simple manifolds. 



The 1^ x Zj 1 ' model on the 3-torua 



The Walker- Wang model based on the Zg 2 - 1 x l!^ 1 cate- 
gory has 2 3 degenerate ground states on the 3-torus, just 
like the 3D toric code. We will not go through the argu- 
ments for this result thoroughly, as it strongly resembles 



those presented for the toric code in Sec. IIIA1 but we 
emphasize that the ground state degeneracy is associated 
with the presence of deconfined excitations. The string 
operator creating a pair of fermions at either end of a 
curve C is: 



n 5? 

crossed edges 



(40) 



where the second product is over edges which are over- 
or under-crossed by C. Here ctJ flips the 2.2 label between 
and 1 on an edge j in state 2.2), while ct| = (— l)* 2 . 

Since Wy commutes with the Hamiltonian, as before, one 
can toggle between the different ground states on the 3- 
torus using the Z 2 operators Wy(C), where C is a closed 
path winding around one of the non-contractible cycles of 
the torus. Furthermore, one can distinguish between the 
2 3 different ground states using the new parity operators: 



p nx = n 



x,y,z 



(41) 



which take values ±1 depending on whether an even or 
odd number of / loops wind around the n-cycle of the 
torus. Here n± is a set of edges parallel to the n-direction, 
whose centers all lie on a plane perpendicular to the re- 



direction. This generalizes the construction in Eq. (111. 
Having understood the Walker- Wang model on a mani- 
fold without boundary, we now investigate its properties 
on a manifold with boundary. 



3. 



The Zg 2 ' x Z^ model on the solid donut 



(2) i 

the bulk, the Z 3 x model is almost identical to the 3 types of string, we can use the reasoning in Sec. |VI C 2 



In the case of the solid donut, the ground state degen- 
eracy is 2 x 3 = 6. Heuristically, this result arises for the 
following reasons. The deconfined fermion can have a 
parity ±1 around the single non-contractible direction in 
the solid donut (this gives the factor of 2). To see where 
the factor of 3 comes from, let us fix the fermion parity. 
This is like considering a Walker- Wang model based on 
the category l}^ . As this category is modular^, and has 
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to show that it has a degeneracy of 3 on the solid donut 
associated with the surface deconfined anyons. There- 
fore, together with the two values of fermion parity, the 
ground state degeneracy is 6. 

Unlike 3DSem, then, the topological order we have 
found on the solid donut does not correspond exactly 
to that of a quantum Hall state on the boundary of a 
confining bulk. The discrepancy is due to the fact that 
we have had to explicitly keep track of the charge e quasi- 
particle in our category, due to its fermionic statistics - 
which effectively doubled the number of ground states. 
The appearance of this deconfined fermionic vertex ex- 
citation is not unexpected: if this state does have the 
topological order of some system of fermions, then the 
fermions should be able to exist in the bulk where they 
would be deconfined in general. 

The extra ground-state degeneracy due to these de- 
confined fermions is, however, somewhat unexpected. It 
arises because the fermions in the Walker- Wang model 
are non-dynamical: The ground state sectors of our 
model can be split into two 'halves', distinguished from 
each other by the number of fermionic strings (mod 2) 
running through the donut. In order to lift this degen- 
eracy we add a term resembling a bulk fermion kinetic 
term t Yltij) ( c l c j + c I c j + ^- c (which also serves to cre- 
ate and destroy pairs of charge e fermions) . In the lattice 
model, the analogue of this term takes the form 



In the case of the fermionic vertex defect, the line-like 
object is not energetic in the lattice model, but is mea- 
surable insofar as it keeps track of any self-twisting in the 
trajectory of the vertex defect. 

The k — 3 Chern-Simons theory certainly captures sur- 
face topological order, but apparently fails to account for 
the additional ground state degeneracy due to the decon- 
fined fermion in the bulk. Furthermore, it is not obvious 
how to use the field theory to create the vortex line de- 
fined in Eq. ( 39 ), or to reproduce the relative statistic be- 



tween the vortex line and the fermionic quasi-particle. To 
capture these additional properties, one might introduce 
another two-component quantum field c^„ and add it to 
the bF + bb action using a term -^e pi,pa c^ u d p A a . The 
world-sheets of the vortex lines (Eq. (39)) correspond to 



(42) 



sources coupled to the new field c via -^c^T^. 

While we have seen that bF + bb theory at k = 3 does 
not describe the lattice model, there is a deep relation 
between the field theory and the category iSp x Z^. 
Chern-Simons theories have a well known relation to chi- 
ral conformal algebras^. Many chiral conformal alge- 
bras, on the other hand, admit a description in terms of 
categories. In the case of k = 3 (more generally odd k) 
Abelian Chern-Simons theory, the co rrespo nding confor- 
mal algebra is in fact a superalgebr a 26 ! 35 !, and the cat- 
egory corresponding to the algebra is Zg 2 -* x iSp ■ The 
Z2-grading corresponds to the fact that Wilson lines car- 
rying label 3 are fermionic under self-twisting; it is pre- 
cisely this Z2-grading which leads to the deconfined bulk 
fermion in our lattice models. 



where ctJ flips the %2 label between and 1 on an edge 
j. This strongly mixes states (ii,0) and such 
that at low energies there remain only three string types, 
which can be labelled by their value of Similarly, 
pairs of ground states differing only by a fermionic flux 
are strongly split, reducing the ground state degeneracy 
from 6 to 3. 



4- The effective field theoretic description 

Given that the surface theory of the Zg 2 - 1 x Z^ lat- 
tice model looks like a v = 1/3 quantum Hall effect, one 
might expect that its topological features are described 
by a bF + bb theory with k = 3 (Eq. J27|). In fact, 
the field theory appears to capture some, but not all, of 
the properties of the lattice model. In the field theoretic 
language, the world-lines of anyons on the surface of the 
sample correspond to Wilson lines carrying (Z 6 ) labels in 
{1, 2, 4, 5}, while the fermions carry label 3; the field the- 
ory reproduces the correct statistics for these particles. 

When we try to source point particles in the bulk, wc 
again encounter the constraint in Eq. ( [29] ): this con- 
straint implies that the particles necessarily lie at the 
ends of line-like objects. This is again consistent with 
the fact that particles labelled {1, 2, 4, 5} are at the ends 
of lines of plaquette defects defects when in the bulk. 



5. More general surface quantum hall states 

Having seen how two specific examples give rise to sur- 
faces with topological order resembling fractional Hall 
states, we briefly summarize how these results general- 
ize. The Walker- Wang model based on category Z^. 1 ^ 
for even k (listed in [26]) gives a fully confined bulk, and 
a surface with the topological order of a v — 1/k bosonic 
Laughlin state. This topological order is captured by an 
abelian bF + bb theory at even level k, written down in 



Eq. (27) 



On the other hand, to obtain the surface topological or- 
der resembling the hierarchical fermionic Laughlin states 
at v = n/m (m odd, n < m and n coprime to m), we 
can construct a Walker- Wang model based on category 
Zm' 1 ' xZj 1 '* (where p is odd and np = 1 mod m, see |26j). 
As in the case of v = 1/3, all of these models have their 
excitations confined in the bulk except for a single decon- 
fined fermionic vertex defect. The lattice models have 
a description in terms of an abelian bF + bb similar to 
that in Sec. |VII B 4| but generalized using the if-matrix 
formalisrrPST 

It is natural to ask whether this correspondence gen- 
eralizes to more exotic Hall states. Using the 12-particle 
category in [55] called 'Moore-Read^ we can produce a 



29 



Walker- Wang model with surface Moore-Read topologi- 
cal order, and a deconfined fermionic vertex in the bulk 
(labelled "0(2,4) m the reference). On the solid donut, this 
model has a degeneracy of 12, but adding the Moore- 
Read analogue of the kinetic term in Eq. (42) lifts the 
degeneracy to 6. 



VIII. CONCLUSION 

In this work, we have discussed two distinct families 
of topological phases in 3D, exemplified by the 3D Toric 
code and 3DSem respectively. The first family (which has 
been studied in a variety of contexts 3 -^) has a ground- 
state degeneracy that depends on the topology of 3D 
space: for example, the ground state is non-degenerate on 
T 3 . The excitations of these models can be grouped into 
point-like 'charges', which are deconfined everywhere, 
and vortex loops; the interesting mutual statistics occurs 
between point charges and vortex loops. These features 
(topological ground state degeneracy, and mutual any- 
onic statistics between point defects and vortex lines) are 
the most natural extension of the notion of topological 
order to 3 dimensions. 

The second family, however, suggests a different notion 
of 3D topological order: it has a unique ground state on 
any closed (orientable) 3-manifold, and all excitations in 
the bulk are confined. Nonetheless on a manifold with 
boundary these models do have degenerate ground states, 
whose number depends on the topology of the boundary. 
More than this, their surfaces admit deconfined anyonic 
excitations. For 3DSem these have (abelian) semionic 
statistics, but other models exist for which these surface 
excitations are non- abelian anyons. Hence these are 3D 
models whose topological characteristics are largely re- 
stricted to their boundaries - but at these boundaries, 
we recover the full richness of possible topological orders 
of 2D systems. 

The interesting low-energy physics of the confined WW 
models thus occurs almost entirely at their boundaries. 
It is well-established that in some systems, such as topo- 
logical insulators 39 , the low-energy physics at the bound- 
ary can be used to classify distinct bulk phases of matter 
(in that case, due to the presence of symmetry-protected 
gapless surface states). We have established that the in- 
teresting surface physics of the Walker- Wang models does 
not play quite so strong a role: since 2D systems with the 
same topological order exist, these surface states do not 
uniquely identify a bulk phase of matter, and indeed can 
be eliminated by adding a 2D quantum Hall layer to the 
surface. 

However, within the realm of "fixed-point" lattice 
models (whose correlation length is smaller than the lat- 
tice constant, and thus effectively 0) the deconfined chi- 
ral surface anyon states cannot be eliminated. Further, 
the existence of chiral anyons at all in an exactly solv- 
able model is notable: we believe that no 2D fixed-point 
Hamiltonian can emulate these^. Thus within the sphere 



of fixed-point Hamiltonians, one can make a strong case 
that the topological order of the surface does uniquely 
identify the bulk. A physical understanding of the mean- 
ing of this bulk-boundary correspondence, which here is 
a feature of the fixed-point Hamiltonian rather than the 
phase, remains an interesting open question. 

Finally, we might ask what kinds of real physical sys- 
tems could be expected to share the long-wavelength 
characteristics of the Walker- Wang models. One pos- 
sibility is that a confined phase of fractional topological 
insulators, where fractionalized excitations are confined 
in the bulk, might in some circumstances be expected to 
support deconfined fractionalized excitations on its sur- 
face. (One of us 4 ^ '■will discuss this possibility in more de- 
tail in a future work.) However, much remains to be un- 
derstood on the subject of realising such phases in more 
physically motivated systems. 

In summary, we have shown that the Walker- Wang 
Hamiltonians present an interesting playground in which 
to examine possible 3D topological states of matter, in a 
context where concrete calculations can be carried out. 
This raises the interesting possibility that some 3D sys- 
tems have a topological order that is purely 2 dimen- 
sional, while others (such as the fermionic models dis- 
cussed here in Sec. VII) are topologically ordered in the 



bulk, but admit much richer anyon models at their sur- 
faces. 
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Appendix A: MTC's are non-degenerate on 
manifolds without boundary 

In this section we show that all confined Walker- Wang 
models have a single ground state on the 3-torus, al- 
though this method of proof generalizes to other simple 
manifolds without boundary. This can also be viewed as 
a proof that there are no deconfined excitations in the 
three dimensional bulk. In Sec. IVI Al we reasoned that 
all confined Walker- Wang models are based on modular 
tensor categories. One of the key features of modular ten- 
sor categories is the 'Flux killing property', described in 
Fig. |23^c): A loop of the string- type w (used to define 
plaquette operators) projects anything it encloses onto 
zero flux. In the following proof we will make heavy use 
this fact, as well as the handle-slide property Fig. [23](a) 
of ujq strings (which holds in all unitary braided fusion 
categories). 
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This proof is structured as follows. In order to show 
there is just one ground state, we need only show that 
V | tp) oc V | 0) holds for any string-net configuration, 
where | 0) is the configuration with no edges colored in 
(meaning all edge bonds have label i = 0) and V is the 
ground state projector. We do this in three parts. First, 
we show that 



v i v) 



E 



(Al) 



where | ip X y)i are configurations with edges colored only 
in the xy plane. We will then apply the same procedure 
to show that, for any of the configurations lying in the 
xy plane, 



V I Ipxy) = / ] bjP 1 VVaxis)i 



(A2) 



where | Vy-axis) are configurations with edges colored 
only along the y-axis. Applying the procedure one last 
time we will show that, for any of the configurations along 
the y-direction 



V | y-axis) = cP | 0) 



Hence, feeding this all into ( Al ), we get 



V | V) « V | 0) 



(A3) 



(A4) 



for all | tp), as required. Let us first prove (Al). Con- 
sider the sequence of diagrams Fig. [27j which show a 
section of the xy plane on the 3-torus. Fig. 27 [a) rep- 



resents some string-net configuration acted on by the 
ground state projector V \ tp) . We remind the reader that 
the ground state projector V = J\ v V v Y\ p V v is a just a 
product of all vertex and plaquette projectors, and that 
a plaquette projector V p can be represented by drawing a 
Wo string in plaquette p. In (b) we act on this state with 
wo loops on the vertical edges sticking up and out of the 
xy plane. We call the operator IC Z , and so we have formed 
a state fC z \ tp). We then act with a set Si of plaquette 
projectors just above the xy plane, each lying in a xz 
plane, to form the state IlpeSi ^p^-^P | tp) in ( c )- Using 
the plaquette projectors in the original expression V \ ip), 
we can handle-slide the plaquette projectors in IlpeSi 
to form the figure in (d), which is another graphical rep- 
resentation of ripeSi 'PpK-z'P I ip)- We can then use the 
flux killing property of wo strings (in MTC's) to show 
that this state is diagrammatically represented by (e). 
We are now near the end of the proof. We act on 

(e) with a set of plaquette projectors just above the 
xy plane to form the state n p eSius 2 ^p^* I V') shown in 

(f) ; we color the new plaquette strings red for clarity. We 
can again use the handle slide properties of the plaquette 
projectors in V | tp) to manipulate the red strings into 
the equivalent form shown in (g). Using the flux killing 
properties of the red strings, we can close up the black 



plaquette strings to form (h). But (h) simply represents 
the state Ilpe™ 'PpP | tp) which is equal to V \ tp). Hence, 
we have proved the following equation: 



II V p K z V\i>)=V\i>) (A5) 

PGS1US2 

Acting on this equation with V, we get the following 
equation: 



vk z v | v) = v | v) 



(A6) 



This equation is the main result of this section because 
the state VK, Z V \ tp) has an interesting property: The op- 
erator JC Z kills all the flux on the vertical legs just above 
the xy plane. This implies that the string-net configu- 
rations in IC Z V | tp) are contractible, at least around the 
z-cycle of the torus. So the string-net configurations in 
K- Z V | tp) can be deformed to lie in the xy plane. Hence 



V}C z V\tp)=J2aiV\tp, 



xy/i 5 



(A7) 



and using Eq. (A6| we get Eq. (Al). We now repeat 



the process. Define an operator K, x by putting w loops 
around edges leaving the yz plane and note the result 
Eq. ( A6 1 also holds for JC X ■ As the operator IC X kills all 



the flux on the legs just to the right of the yz plane, the 
string-net configurations in K, X V \ tp xy ) are contractible 
around the x-cycle of the torus (as well as the z-cycle). 
So the string-net configurations in JC X V \ tp xy ) can be 
deformed to lie along the y-axis. Hence 



VK, x V\tp xy )=Y J hV\tp v 



(A8) 



and using Eq. ( A6 1 we get the result Eq. ( A2 ) . Now to 



the final stage of the argument. Define an operator JC y 
by putting ujq loops around edges leaving the zx plane 



and note the result Eq. (A6| also holds for JC y . As the 



operator K, y kills all the flux on the legs leaving the zx 
plane, the string-net configurations in Yi y V \ y-axis) are 
contractible around the y-cycle of the torus (as well as 
the z-cycle and x-cycle). So the string- net configurations 
in JC X V | VVaxis) can be deformed to | 0). Hence 



VK X V | i/y. 



cV\0), 



and using Eq. (A6) we get the result Eq. (A3), 
statement Eq. (A4) follows from this. Thus we have 



(A9) 



The 



shown that the flux killing property of MTC's implies 
that Walker- Wang models based on MTC's are non- 
degenerate on the 3-torus. 
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FIG. 27. (Color online) This figure shows the sequence of steps used in proving that Walker- Wang models based on MTC's are 
non-degenerate on the 3-torus. 

Appendix B: All string operators are confined in plaquette violations along their length. This result sup- 

MTC's plements the discussion in the text, where we considered 

a special type of string operator (Wy(C)) and showed 

t ,i • ,i , ( ,1 ..mn, , that it fails to commute with the plaquettes along its 

m this section we prove that for all MIC s, any vertex ^ ^ 1 H 6 

types string operators (obeying mild constraints) produce en ^ 
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Suppose we construct a string operator V(C) with the 
following qualities. It can change the color of edges along 
a path C, and assign a configuration dependent phase to 
edges touching the path C. Furthermore, we assume that 
if the edges on C have label i = 0, V(C) will flip the edges 
to some superposition of states, some of which have all 
the edges on the path C in a state i 7^ 0. In this section 
we will show that such an operator necessarily fails to 
commute with plaquettes along its length. 

Suppose that C pierces the xy plane once. Take a large 
set S of plaquettes in the xy plane with the condition that 
C pierces the set, but the plaquettes n ear the boundary of 



K(dp) 



S are not acted on by V(C) (see Fig. 28). We will show 
that, with the above conditions, V(C) cannot possibly 
commute with some plaquettes in the xy plane. This 
shows that the defects created by V(C) are necessarily 
linearly confined. 

With these above assumptions about V(C), we can say 
that upon acting on the empty ket | 0) (all edges carry 
label i = 0), the string operator flips all the edges along 
C to some superposition of labels s — J2i 7^ 0- With- 
out loss of generality (by subtracting a multiple of the 
identity operator from V(C)) we can assume that so = 0. 
This implies that 



W$°(dS)l[V p V(C)\0) 

pes 



(Bl) 



where Wy°(dS) is a string operator with the special 
handle-slide property, and dS is some path on the bound- 
ary of the set of plaquettes in the xy plane (see Fig. 28). 



This equation follows simply from the fact we can handle- 
slide the perimeter string dS inwards so that it encloses 
the edges on path C which carry label s. As we are in a 
MTC based model and sq = 0, we can use the flux killing 
property of the Fig. [23^b) to get zero in Eq. (Bl ). 

If we assume that ri P eS an( ^ ^(^) commute then we 
are led to a contradiction. Continuing from the previous 
equation: 




FIG. 28. (Color online) This figure represents a string op- 
erator V(C) piercing the xy plane. The black loops are the 
plaquette operators in a set S, pierced by the string operator, 
while Wy°(dp) (blue line) is another string operator (with 
label (jJq) which encircles the perimeter of the set S). In the 
text we show that V(C) fails to commute with some of the 
plaquettes in S, with fairly mild assumption on the form of 
V(C). 



Therefore, 



= V(C) Y[v P \o) 

pes 

= n w c ) 1 °> 

pes 



(B3) 



But this state is just a cluster of plaquette projectors 
acting on a ket with some edges in superposition s. It can 
be shown that such a state is not zero by expanding out 
the plaquette projectors in terms of plaquette string op- 
erators Wy(dp). Consider a string configuration where 
all edges have label except for those on C (which are 
in superposition s). This configuration can be shown 
to occur once in the expansion of Yl peS V p V(C) | 0), 
and with a non-zero coefficient. Hence the expression 
Ilpgs 'PpViC) I 0) cannot be zero. Thus, the assumption 
that YipeS^P an( ^ y(^) commute leads to a contradic- 
tion, implying V(C) is linearly confined. 



= W^(dS)l[V p V(C) |0) 



pes 



W^(dS)V(C)]Jv p \0) 



pes 



= V(C)W^(dS)l[V p \0) 

pes 



(B2) 



To get to the second line, we assumed that the pla- 
quettes commute with V(C). In the third line, we used 



the fact that 



V(C),W^°(dS) 



because they do not 



act on the same edges. Now, Wy° (dS) ]\ p£S V p \ 0) = 
V p I 0) because we can handle-slide the string over the 
plaquette operators to get a wq loop enclosing no flux. 



Appendix C: Explicit forms for operators in DSem 
model 



In our main treatment of the DSem model, we defined 
the plaquette operators only in the absence of vertex de- 
fects. Furthermore, we never explicitly showed the form 
of string operators near their end-points. In this ap- 
pendix, we address some of these issues. In the main 
text we defined the plaquette operators as 



Bp-di^ n * (i ^ )/2 

iedp jes(p) 



(CI) 



The difficulty with Eq. (CI) is that as operators, B pi 
and B P2 do not commute when acting on neighboring 
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FIG. 29. The figure shows an adjacent pair of plaquettes (pi 
and P2), used to show that adjacent plaquette operators in 
the DSem model (as defined in Eq. Q) do not commute if 
vertex defects are present. 



plaquettes. Specifically, on the 5 edges on which both 
operators act (see Fig. 29 1 we have 



Bp x B P2 — 



[a 1 CT3 <7r 



B P2 B pi 



(C2) 




(A) 




FIG. 30. (a) Shows a labelled plaquette, used to define the 
complete plaquette operator, (b) Shows a graphical represen- 
tation of the string operator defined in the text. 



We have defined Q v = {I-B v )/2, and A„ = J2 e es(v) n ^ 
where e runs over the edges connected to vertex v e.g. 
Ai = riQi + nil' + ^12- This plaquette term represents 
a rather complicated twelve spin interaction term, al- 
though we remind the reader that this reduces to the 
somewhat simpler form (the first line of Eq. (C6)) in the 
absence of vertex defects. 



where 



B Pl ,B P3 ] = -iafa$a%a% sin 6*1245 



?1245 = ^ (02 + °\ ~ °1 ~ a l) 



In the ground state, where 

n 



B v = [[ a! I 

iedv 



(C3) 
(C4) 

(C5) 



6*1245 is a multiple of 7r, and the two products are equal. 
Thus the plaquette operators (CI) commute only when 



acting on states without vertex violations. In states with 
vertex violations, 6*1245 is an even or odd multiple of 7r/2, 
depending on whether there is an even or odd number of 
vertex violations on the two vertices (123), (345). 

For the following definition of B p it is more convenient 
to label the vertices of p rather than edges (Fig. 30 ^a)), 
and denote edges by the pair of vertices connecting them. 
Edges are now labelled by the pair of vertices they con- 
nect: 



B p = II o?i+i(-l) n *- 1,(1-n " +0 
j-(l-Qi)"ii'Ai^(l-Q2)n 22 /A 2 

j(l-Q3)(l-n 3 3')("23-n34) j-(l-Q4)«44' A 4 
j(l-Q5)"55' A 5 j(l-Q6)"66'( n 56-nei) 



(C6) 



1. String operators, and their endpoints. 

In this section we present a formalism for introducing 
pairs of vertex defects into the DSem model. We will 
write down the form of an open vertex-type string op- 
erator, near and away from its endpoints. We use the 
convention that positive chirality particles are associated 
with an over-crossing type string, and negative chiral- 
ity under-crossing type string. In the formalism used, 
positive chirality particles will also be associated with a 
plaquette defect; in fact the most natural formalism re- 
quires that we introduce an additional label at vertices, 
but we opt to use another convention where one chirality 
of string involves a plaquette and vertex defect at the 
end-point. 

For Fig. [30{h), we make the negative chirality string 
operator take the following form if acting on a state with- 
out vertex violations: 



ik n (- 

iec <jfc)ec 



1) 



nj(l-n fc ) 



n • 

1) at R vertex 



•JA(PB ■ 



(C7) 

Here the 4>a,4>b factors are associated with the precise 
state of the string end-points. The phases are defined 

i (l^n Ar )n Al whne ( j )B = (_^(l-n Br )n Bl ^ where 

> Bi, 



9A = 

Ai,A 



The 



B r are the edges marked in Fig. 30 'b). 
reader should be warned that the exact forms of these 
phases will depend on which direction the string ap- 
proaches the vertex, as well as the convention adopted 
for incorporating vertex defects. 

To produce a defect of opposite chirality, we conjugate 
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the operator, except (pA,<f>B parts. The resulting oper- 1 2 



ator is of positive chirality, and will anti-commute with 
the plaquette labelled p at its endpoint, but no other 
plaquettes. 




Appendix D: Topological entanglement entropy in 
3DSem 



Here we show that unlike its 2D counterpart, the 3D 
semion model has no topological entanglement entropy. 

As described in Ref. 2U there are two sensible pre- 
scriptions for defining topological entanglement entropy 
in 3D. We will use one of these for our calculation, but 
the result is independent of the choice of prescription. 
Both prescriptions operate on the same principle, gener- 
alizing the approach of Levin and WerP^for 2D systems: 
we add together the entanglement entropies of several 
possible partitions of a system into subsystems A and 
B, with coefficients such that the net boundary and cor- 
ner terms all cancel. Fig. [31] shows the combination of 
partitions that we will use. We will always cut the two 
subsystems along the middle of a set of edges, so that the 
two subsystems must share a set of edge labels along the 
boundary. 

It is useful to briefly recall the source of the topological 
entanglement entropy in the Toric code. The Toric code 
ground state is a superposition of all configurations of 
closed loops, with a relative amplitude of 1. Thus, for 
any partition of the system into two subsystems A and 
B, we may write 



i*o> = £«d*A e) >i*2 e) } 



(Dl) 



where i c denotes a particular choice of edge labels on the 
boundary between A and B, for which the total number 

of occupied edges crossing the boundary is even. ) 

and I'F^) are themselves superpositions over many dif- 
ferent loop configurations, with the configuration i c of 
spins on the edge. |^o) factorizes according to Eq. ( |D1[ ) 
because the relative coefficient of all elements in these 
superposition is always 1. 

Because configurations in which the edge labels are 
different are orthogonal, it is easy to compute the reduced 
density matrix: 



PA = Tr B 



E^(^)i^i 2 i^ c) )M o) i 



where -/Vb(z c ) is the number of configurations in B with 

these boundary conditions. Since \^a^) x &cd, 

Pa is diagonal, and we may read off the entanglement 
entropy: 



>AB 



= y / N B (i c )N A (ic)\ai c \ 2 log [N B (i c )N A (i c )\i 




FIG. 31. The combination of partitions used to calculate the 
entanglement entropy. The figure shows the boundary of the 
regions A and B. 



subject to 



,N B (i c )N A (ic)\ai 



(D3) 



In practise NA{i c )^B{i c ) will be independent of the 
particular boundary configuration, as will |aij 2 , so that 
NA{ic)NB(i c )\oti a \ 2 = 1/M, where Mi is the number of 
possible boundary conditions. In a loop gas, if the total 
number of edges on the boundary between A and B is 
tu, then 



Mi = T 



-N 



(o) 



(D4) 



where is the number of connected components of 

the boundary of region A. The entagnlement entropy is 
thus : 

S = ]T jjr log M t = (n, - < } ) log 2 (D5) 



The topological entanglement entropy is given by the 
combination of subdivisions shown in Fig. |31| In the first 
three, A consists of a single component and N A — 1; in 
the last term — 2. Since the regions are chosen such 



that 



,(2) 



,0) 



0, we obtain: 



'Top 



(n 



(i) 



(n 



(2) 



Of ) - 1) + (n\ 4) - 2) log 2 = log 2 



For the doubled semion model, if A is simply con- 
nected then we may use exactly the same reasoning as 
for the Toric code to calculate the entanglement entropy 
(whether or not A has multiple boundary components). 
That is, in this case the ground sta te w ave function can 



be decomposed according to Eq. (Dl), with 



(D2) superposition of all loop configurations in A with the 
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FIG. 32. The scenario if A is not simply connected. If the 
total number of down spins crossing the cut through the donut 
is odd, then there is a loop that is not contractible within 
region A. Similarly, if the number of down spins through the 
center of the donut is odd, then there is a loop that is not 
contractible within region B. When both of these loops are 
present, the wave-function has an extra phase —1 due to their 
linking. 



boundary spin configuration i c , and similarly for \^f B ). 
Explicitly, when A is simply connected every loop con- 
figuration in A with a fixed boundary i c can be obtained 
from every other such configuration by acting with some 
number of plaquette projectors on plaquettes that are 
entirely inside A. Unlike in the Toric code, there are 
non-trivial relative phases between these different loop 
configurations; however, these phases are entirely dic- 
tated by the action of the plaquette projectors within A, 
and are independent of the loop configuration in B, such 
that the wave function still factorizes. The rest of the 
calculation is identical to that for the Toric code; hence 
if A is simply connected, the entropy is again given by 
Eq. JD5]). 



If A is not simply connected, in general the factoriza- 
tion (Dl) fails because the doubled semion ground state 



contains phases that are sensitive to the linking number 
of loops in A with loops in B. The case of interest is 
where A is topologically a donut. Let us separate the 
wave function |5'^ c ' 1 ) according to: 



I* 



(»c)\ 



I* 



(0\ 



(D6) 



ber of down spins through the hole in the middle of the 
donut is even or odd (i.e., according to whether or not 
there is a loop in B that is not contractible within B): 



I* 



(ic)\ 



1 

7^ 



B )o 



I* 



(ie)\ 



(D7) 



The ground state of 3DSem is given by 



E 



a. 



I* 



( i c)\ 



I* 



(ic)\ 



I* 



I* 



B 1° 



I* 



(i c ). 



(ic)\ 



(D8) 



The extra — sign in front of the last term (in addi- 
tion to any relative phases present in the definitions of 

|*^ c ^)e,o, |^B^)e,o) arises because when both odd com- 
ponents are taken, the non-contractible loops in A and B 
are linked, which produces an extra — sign in the wave 
function. 



Because of this extra sign, Eq. (D8| cannot be ex- 



pressed in the form Eq. (Dl I; rather, we have 



l*0)=E^{(l^ 3 )e + l*^ ) )o)|^ ) 
l^o) 1^)0} 



so that 



{(' 



(ic)\ 



-N B (i c ) 



(D9) 



(D10) 



^>o) ((^le + M^lo) 
+ (l*S° ) >.-|*2' ) >o) ((^le-^lo)} 

Here we have used the fact that the number of loop con- 
figurations in the even and odd sectors of | v t / ^ c ' > ) are 
equal. As this also holds for |^/^ '), we have: 



I* 



(»c)\ 



(Dll) 



so that 



Tr PA = J2 W\ 2 N A {i c )N B {i c ) = 1 (D12) 



Here e and o refer to whether the total number of down 
spins measure d a long a slice through the donut is even 
or odd (Fig. 



32). |1^ c) ) c (and similarly |*^ c; )o) is a 



(«e)\ 



superposition of many different loop configurations with 
different relative phases, which are related by the action 
of some number of plaquette projectors acting entirely 
inside A. To get from a configuration in the even set 
to one in the odd set, however, one must act with pla- 



(ic)\ 



and 



quette projectors in region B as well. Thus \^ A 

|\l/^)e need not appear in the ground state with the 
same phases, for a given loop configuration in B. 

To see the implications of this, let us separate the loop 
configurations in region B according to whether the num- 



Since the coefficients \atij\ 2 N A(i c )N b{i c ) ar e independent 
of the choice of boundary spin configuration, this fixes 



a ic \ 2 N A (i c )N B (i c ) 



1 

Wi 



(D13) 



where ni is the total number of edges crossing the bound- 
ary between A and B. This gives the entanglement en- 
tropy: 



Sab = -J2\ 



2 N A (i c )N B {i c ) log 



1 



2 N A {i c )N B {i c ) 



■ Mijjr \og2Ni = Hi log 2 



3G 



31 are simply 



for the solid donut. 

Since the other three regions in Fig. 
connected, their entanglement entropy is identical to that 
of the toric code. The topological entanglement entropy 
of the 3D semion model is thus: 



STop — 



K W ) - («r - 1) - W 1 1) + (nf> 2) log2 



,(2) 



,(3) 



(4) 







exactly as one might have expected for a system that has 
no bulk topological order. 



Appendix E: Surface plaquette operators of the 3D 
semion model with rough boundary 

The choices of boundary conditions in 3D are similar 
to those in 2D, which have been discussed at length in 
Refs. [T8l and fT9l In the main text we always assumed a 
"smooth" boundary. Here we consider a "rough" bound- 
ary where "dangling" edges are sticking out of the sur- 
face. Including these these dangling edges on the bound- 
ary allows a source of a colored loops at the boundary 
with no energy cost. 

For the toric cod d 18 * 19 -* one can create "partial- 
plaquette" operators which flip the value of two neigh- 
boring dangling edges as well as the bulk edge connect- 
ing them. These partial-plaquette operators commute 
with the Hamiltonian and have the effect of allowing the 
sources to move freely along the surface. 

The situation with the 3DSem model is a bit differ- 
ent. We will similarly be able to construct such partial- 
plaquette operators which commute with the bulk Hamil- 
tonian, however, they will not commute with each other. 
This is to be expected since the sources should have 
semionic statistics with respect to each other. 

The form of such a partial plaquette is show in Fig. 33 
The plaquette operator acts by flipping a pair of neigh- 
boring dangling edges (labeled 1 and 13 in the figure), 
together with the bulk edges connecting them (3, 5, 7, 9, 
and 11). In order that this operator commute with the 
bulk plaqucttes that also act on these bulk edges, we 
must include the usual phases depending on the edges 
bordering those that are flipped (2,4,6,8,10, and 12 in 
the figure), together with an extra phase depending on 
the value of some of the flipped edges before flipping. (In 
the bulk plaquette operator (14 1, these are the red and 
blue edges 



For the surface plaquette drawn here, only 
one red edge (3) and one blue edge (1) are included). 
This gives 

where 



(E2) 



However, there is a problem here: B p0 is not Hermitian 
if m — ni3 is odd. To compensate for this, we must add 




A 



FIG. 33. (Color online) A surface plaquette. The surface 
plaquette operator described in the text flips the thickened 
edges (1, 3, 5, 7, 9, and 13). It has a phase that depends on rii 
on he edges bordering those that are flipped (2, 4, 6, 8, 10, and 
12, shown here in green), as well as an extra phase depending 
on the occupancy of link 3 (red) and link 1 (blue). In this 
figure the dangling edges are 1 and 13. 



an extra phase factor (— l)" 1 ™ 13 . The final form of the 
surface plaquette operator shown is: 



psurf _ x x x x x x ■n 2 +n 4 +n e +n8+n 10 -n 3 -ni ( i \r 



More generally, we have 



(E3) 



b s x = ( n a ^ n *™^ Ej r ° d ™ 3_E " wuo " j dans " ni 

iedp j£s(p) 

(E4) 

with red and blue edges defined differently for the differ- 
ent plaquette orientations, as in the bulk. 

In Sec. IIV B 21 we claimed that we could add a kinetic 
term to the vertex defects on the surface of the manifold. 
This is done by adding dangling edges and including these 
partial plaquette operators. 



Appendix F: Trace methods 

In this section we give a more standard derivation 
of the ground state degeneracies of the toric code and 
3DScm models: We will prove that the 3D toric code 
has ground state degeneracy 2 bl on a manifold without 
boundary (and b\ independent non-contractible cycles), 
whist the 3D semion model has a single ground state. 
For example the 3D toric code and semion model on the 
3-torus have ground state degeneracies 2 3 and 1 respec- 
tively. In Sec. |F 1| we introduce a standard method used 
to calculate ground state degeneracy and apply it to the 
2D models. Then, in Sec. |F 2| we generalize the method 
to the 3DSem model. 



1. 2D ground state degeneracies via trace method 

In this section we algebraically deduce the ground state 
degeneracy of the toric-code and DSem model embedded 
on some 2D surface without boundary; we have already 



37 



discussed the ground state degeneracy using a loop gas 
picture in Sec. [TTJ but we opt here for a more concrete 
method. This section should be viewed as preparation 
for the work in 3D in the next section, where we use a 
similar method. We define a projector onto the ground 
state 



v = n v « n v p 



(Fl) 



peP 



where V is the set of all vertices, and P is the set 
of all plaquettes. In the toric code V v = (1 + B p )/2 
where B p is defined in Eq. ([!]). In the DSem model V p = 
(1 — B v )/ 2 where B p is defined in Eq. Q. In both models 
V v = (1 + B V ) is the same. To calculate the ground state 
degeneracy Tr (V) we write the trace as 



Tr(V) 



F x 2 x J2 Tr 

SvQV 



n & 

.v£S v 



(F5) 



We use a similar line of reasoning to show that only two 
terms survive the remaining sum. Note that 1C V is -1 if 
vertex v has an odd number of 1-strings entering it, and - 
1 otherwise. Therefore the product H^eSv ac ^ s * n ^ e 
following way. Edges that do not belong to a vertex in 
Sy are unchanged, and those that belong to two vertices 
in Sy are also unchanged because (a x ) = 1. However, 
those boundary edges belonging to just one vertex in Sy 
are acted on with <r x . The fact that a x is traceless im- 
plies that Tr [Ilues K-u] = unless Jl^es^ ^ v ^ as no 
boundary edges i.e. unless Sy = ° r V. In both cases 



n 



ves v 



JCv = 1, therefore 



Tr (V) = 2 



-V-F 



l[(l + IC v )l[(l + IC p ) 



p£P 



= 2 



-V-F 



E E Tr 

SvQV SpCP 



Y\_ fcv JJ fc t 



vev 



peP 



(F2) 



(F3) 



where K v — 2V V — 1 and K, p = 2V P — 1. Here Sy and 
S p are subsets of the set of all vertices and plaquettes re- 
spectively. For the toric code and DSem model, all of the 
terms in the sum are zero except for the four correspond- 
ing to (Sv,S P ) e {(0, 0) , (V, 0) , (0, P) , (V, P)}. For each 
of these 2 2 = 4 choices, n^eSv ^« IlpeSp = anc ^ so 



Tr{V) 



-,2-V-F 



Tr[l] 



2 2"X 



Tr (V) = 2 2 - v - F Tt [1] 
= 2 2 ~ x 



(F6) 



where x 1S the Euler characteristic of the surface that 
the lattice tiles, and so the ground state degeneracy of 
the toric code is is 2 2g where the integer g is the genus of 
the surface. Indeed the answer is the same in the case of 
DSem because the IC V are the same and the K, p = —B p 
still flip edges; the only non-trivial fact to check is that 

UpeP^p = L To act with UpeP^p = (- l ) F Y[peP B P 
we can first fuse the B p operators on each edge (which 
gives a factor of —1 for each edge). All that remain are 
the closed 1-loops on each vertex, which give a factor of 



-1 for each vertex. In total, Yiptp^p 



■1) 



V-E+F 



(-1) 



2a 



1. Hence we have proved the statement in 



(F4) Eq. fF6| 



where X = V" — E + F is the Euler characteristic of the 
surface tiled by the lattice. Therefore the ground state 
degeneracy of the 2D toric code, and DSem is 2 2 ~ x = 2 9 
on a surface of genus g. 

To see how the above argument works in detail for the 
toric code note that JC P = B p acts on each edge of pla- 
quette p with a a x i.e. it flips the spin on each edge 
of p. Therefore the product IlpeSp ^p acts m tne 
lowing way. Edges that do not belong to a plaquette 
in Sp are unchanged, and those that belong to two pla- 
quettes are also left unchanged because they are flipped 
twice. However, those edges belonging to just one pla- 
quette in Sp are flipped once by IlpeSp ! we wm 
call these boundary edges. If the operator has boundary 
edges then for any state | ip) in our a z = ±1 spin basis 
I UveSv Upes P I V>) = 0, because U peSp IC P 
will flip at least one edge in | ip) to an orthogonal state, 
and Ilties , cannot undo this. Therefore we need 
only worry terms for which JlpeSp nas no boundary 
edges, which occurs when Sp = or P. In both cases 
Ilpe s P ^p = 1 > therefore 



2. 3D ground state degeneracies via trace method 

Having calculated the ground state degeneracies of the 
2D models, we now apply similar methods to the 3D 
toric code and 3DSem on manifolds without boundary. 
We will confirm the results presented in Sec. |III AT] The 
3D toric code has degeneracy 2 bl where b\ (the first Betti 
number) is the number of independent non-contractible 
cycles in the closed 3-manifold, while the semion model 
has a single ground state. 

Following our work in 2D Sec. |F 1| we cast the ground 
state degeneracy in the form 



Tr (P) = 2 



-V-F 



Y\_ fc-v JJ ic p 

vev peP 



(F7) 

where JC V — 2V V — 1 and JC P — 2V P — 1. We evaluate 
terms in the above sum, first for the toric code and then 
for the semion model. 
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a. 3D toric code: }C P = B p flips the spin on each 
edge of plaquette p. Therefore the product IlpeSp ^p 
acts in the following way. Edges that do not belong to 
a plaquette in Sp are unchanged, and those that belong 
to an even number of plaquettes are also left unchanged 
because they are flipped twice. However, those edges 
belonging to an odd number of plaquettes in Sp are on 
net flipped once by the operator Yl P e Sp^p> we again 
call these boundary edges. If the operator has boundary 
edges then for any state | ip) in our <j z = ±1 spin basis 



(vi n ^ n £pi^>=o, 



(F8) 



v eSv 



peSp 



because IlpeSp ^p ^ip a ^ least one edge in | tp) to 
an orthogonal state, and Ilues^ ^ v cannot undo this. 
Therefore (F8) is zero, except possibly when Sp has no 



boundary edges. It can be shown (see below) that there 
are 2 bl+c ~ 1 sets Sp with no boundary edges, where 
C is the number of cubes on the lattice. Furthermore 
Yisp fcp ~ 1 f° r Sp with no boundary edges. This leaves 
us with a ground state degeneracy of 



start again with ( F7 1 recasting it as 



Tr{V) = 2' 



Tr 



SpCP 



n v * n 



veV 



pes F 



(F10) 



The K p = —Bp operators in the 3DSem flip each edge of 
the plaquette p (albeit with possible phases), so we can 
use the same reasoning as for the toric code to show that 



ripeSp ^p IS traceless unless Sp has no boundary edges. 
We found before that there were b 2 types of Sp without 
boundary edges by remembering that any such Sp can 
be visualized as a set of possibly intersecting surfaces 
without boundary. This has not changed. However we 
will find that only one of these b\ types of surface gives 
non-zero trace: IlpgSj, is traceless unless Sp has no 
boundary edges and can be represented as a contractible 
surface. 

To see how this works we recall the result shown in 



Eq. (18) that P n(Sp )± = Upesp (- R p)> where P * 



n(Sp)± 



IS 



the parity ±1 of the number of loops crossing the non- 
contractible closed surface defined by Sp, which we say 
has normal n(Sp). Therefore 



Tr (V) = 2 



-v-F x 2&1+C 



1 >< E Tr 



s v cv 

2 -V-F x 2 b 1+ C-l x 2 1 + E 
2 bi-X 

2 bl 



n ^ 



(F9) 



where x = C — F + E — V is the Euler characteris- 
tic, which is zero for closed 3-manifolds (Poincare dual- 
ity). To show the second equality we used the fact that 
fives ^ v i s t race l ess unless Sv — or V (see Sec. F 1 1. 
In both cases ELes K, v = 1. 

To justify the factor of 2 bl+c_1 above, note that a set 
Sp without boundary can be visualized as a set of pos- 
sibly intersecting surfaces without boundary. Consider 
first those Sp with surfaces that can be contracted to a 
point. Such an Sp can be said to form the boundary of 
one of two sets of cubes: Sc or its complement C — Sc, 
where C is the set of all cubes. Therefore those Sp whose 
surfaces can be contracted to a point can be put into 1:2 
correspondence with the set of subsets of C, and so there 
are 2 C ~ 1 such Sp. But not all Sp have surfaces that can 
be contracted to a point. In fact there are 2 &2 topologi- 
cal classes of surface where 62 is the second Betti number 
(which for oriented closed three manifolds obeys 62 = b\). 
Each of these classes of surface can be shown to contain 
2 C ~ 1 elements in the same we showed above for surfaces 
that could be contracted to a point. 

b. 3D semion: Having found that the ground state 
is degenerate for the toric code on closed 3-manifolds, 
we now show that the ground state is non-degenerate for 
the semion model on any oriented closed 3-manifold. We 



Tr 



n v » n ^ 

vev peSp 



Tr 



i I P V P, 



(Sp) 



.vev 



(Fll) 



However, consider the invertible operator S = 
Y\iec n a i> wner e Cn is a closed path cutting Sp perpen- 
dicularly. Notice that S commutes with V v but anti- 
commutes with the parity P n (Sp)± so that 



Ti- 



ll V - P r 



(Sp) 



.vev 



Tr 



S 1 Y\_ p vP n {Sp)±S 



vev 



= - Tr 



n 

vev 



PvP n (Sp)j 



(F12) 



Hence, the trace must be zero for Sp. If Sp is con- 
tractible, then an odd number of fluxes pierce it only if it 
encloses vertex violations, but such states are projected 
out by U veV V v . Therefore 



Tr (V) = T F x ^ Tr 

s P es 

= 2- F x \£\ x Tr 



n v - n k v 

vev peSp 



vev 



= 2- F x x 2 1 - v Tt[1] 

= 2~ x 
= 1 



(F13) 



where x is the Euler characteristic (which disappears 
for closed 3-manifolds), and £ is the set of closed and 
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contractible 5 p 's which has size 2° (as reasoned in the case of the toric code). 
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